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Multiple peak aggregations for the Keller-Segel system 
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P\J , Abstract 

^ , In this paper we derive matched asymptotic expansions for a solution of the Keller-Segel 

system in two space dimensions for which the amount of mass aggregation is SttTV, where N = 
1, 2, 3, ... Previously available asymptotics had been computed only for the case in which A^ = 1. 

^^ . 
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^ ; 1 Introduction 

-(— > 
c^ _ 

The goal of this paper is to describe, using matched asymptotics, the asymptotic behavior near 

blow-up points of a class of nonradially symmetric solutions of the following Keller-Segel system. 



^; ut = Au-V (uVv), X € M^, i > 0, (1.1a) 

^ = Au-Fu, xGM^ t > 0, (1.1b) 

Q , The Keller-Segel system, which was introduced in [12], is a classical model of chemotactic aggre- 

^r^ I gallon. In this model u is the density of a biological organism and v is the concentration of a 

^D ■ chemical substance produced by it having chemoattractant properties. It was conjectured in [3] 

^^ and rigorously proven in jTT], in the case of bounded domains, that solutions of (ll.lap - (ll.lbp may 

blow-up in finite time, showing the fact that is usually interpreted as the formation of a high density 
aggregate of cells. 

/\ . The mathematical properties of (|l.la|) - (|l.lb[) have been extensively studied by many authors, 

j^ I One of the most peculiar features of ()l.lap - ()l.lbp is the existence of a critical mass ttiq such 

that for solutions with initial total mass of organism J uq larger than ttiq, blow-up takes place, 
whereas solutions with smaller values of J uq yield global existence of solutions (cf. [H [T3] , for 
bounded domains, [U [15] in the case of M^). It has been already proven that blow-up consists in 
the formation of a Dirac mass in finite time with an amount of mass larger than An in the case 
of Neumann boundary conditions and blow-up taking place at the boundary of the domain, and 
larger than Svr in the case of blow-up taking place at interior points (cf. [17]). The literature about 
the Keller-Segel system is huge and we will not attempt to summarize hear all the existent research 
concerning singularity formation and global existence for (jl.lap . (jl.lbp . Some of the main results 
in this direction can be found in [T| [2t 14} [TT | [Ti l [TS] . 
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In the case of radially symmetric solutions, the asymptotic behavior of solutions of (jl.lap - (]l.lbp 
near blow-up points was obtained in [8] using asymptotic methods, and a rigorous construction of 
such solutions was given in [9]. Actually the paper [8] describes formally the asymptotics of the 
blow-up solutions also in the parabolic-parabolic case in which (ll.lb|) is replaced by a parabolic 
equation. The rigorous construction of the corresponding solutions is given in [10]. The solutions 
constructed in [9] produce the aggregation of a Dirac mass with the mass Svr. On the other hand, 
continuation of solutions after blow-up has been considered using formal arguments in |201 [2T] , and 
rigorous mathematical analysis in [U [13] . 

We will describe in this paper the asymptotics of solutions of (|l.lap , (|l.lbp yielding formation 
of Dirac masses whose amount of mass is SvriV with A^ = 2, 3, 4, ... These solutions will be obtained 
by means of the coalescence at time t = T oi N peaks of mass placed at distances of order 
\/T — t, each of the peaks containing an amount of mass asymptotically close to Svr. The behavior 
of such solutions will be obtained using matched asymptotics. The peaks where most of the mass 
is concentrated near the blow-up time are placed at the vertices of some polygons to be described 
in detail later. 

We summarize the main result of this paper in the following Theorem. We emphasize that the 
results of this paper are obtained at the level of formal asymptotic expansions but not a rigorous 
Theorem in the sense of Mathematical Analysis. 

Theorem 1.1 It is possible to find formal asymptotic expansions for solutions of the Keller-Segel 
system lll.la]) . 1(1. lb\} that blow up at the time t = T at the point x = xq and at each time t < T 
the mass is concentrated around the points Xj (t) , j = 1,2 where: 



Xj (t) = xo + i-iy^^ aVT - t, j = 1, 2 

with a = (2, 0) G M^. More precisely, the formal solutions described by the asymptotics found in this 
paper have the following property. For any ly > arbitrarily small, there exists R > sufficiently 
large such that: 

lim / u {x, t) dx — 8iT < v 

with 

for some a > 0. 

Moreover, the total amount of mass concentrating at the point x = xq as t ^ T~ is 16tt. More 
precisely, for any function 7] (t) such that lim(_^7i- r] (t) /\/T — t = oo and liui^^'j^- r] (t) =0 one 
has: 



lim / u (x, t) dx = IGvr. 



'B,,(t)(a:o) 

Remark 1.2 The argument used in the construction suggests that it would be possible to obtain 
solutions yielding the aggregation of an arbitrary number of multiples of Sir. However, the feasibility 
of such a construction requires to check that a certain elliptic problem, associated to suitable singular 
self-similar solutions of (jl.lap . (jl.lbp (cf. Section \^, satisfy some sign condition that will be 
discussed in detail in Section [7] for the case in which two peaks aggregate. We have checked that 
this sign condition holds in this particular case solving numerically an elliptic equation. Analogous 
sign conditions should be checked for aggregations of multiple peaks, which we have not attempted 



in this paper. Precise asymptotic formulas for the solutions described in Theorem \1.1\ will be given 
in the rest of the paper. In particular, we will derive precise formulas for the width of the regions 
around the points Xi (t) where the mass concentrates. The final profile of the solution at the blow-up 
time will be described in Remark \5J 



The results of this paper are of a local nature. For this reason we just restrict our analysis to 
the case in which the system is solved in the whole M^. Similar results could be derived for the 
Cauchy-Neumann problem in bounded domains with non-flux boundary conditions (cf. Section [7]). 

We finally remark that numerical simulations showing aggregation of several peaks at the time 
of the singularity formation were obtained in J16j . 

2 Notation and preliminaries. 

As indicated in the Introduction we will denote as T the blow-up time. We will use repeatedly in 
the rest of the paper the following self-similar variables: 

w(x,i) = ^-^$(y,T), v{x,t) = W{y,T), (2.1a) 

y = ^=^, r = -log(r-t). (2.1b) 

The system (jl.lap . (|l.lbp becomes in these variables: 

$^ = A$-^ V($VVF)-$, (2.2a) 

= AW + ^. (2.2b) 

It is natural to expect a self-similar behavior for the solutions of (I2.2ap . (|2.2bj) . Self-similar solutions 
of (fTTill . (fLTbjl solve: 

A$ - y^ - V ($VVF) - $ = 0, (2.3a) 

AW + ^ = (2.3b) 

in the variable (j2.1ap . (|2.1bp . The solutions that we construct in this paper approach asymptotically 
as T — )• oo the singular steady states: 

N 

^s = 87rY,Hy-ye) (2.4) 

with the points yi satisfying: 

N 

1-4 Y. ^^ = 0' J = h2,...,N. (2.5) 

The solutions ()2.4p . (12. 5p solve p.3ap . (I2.3bp in the sense that they can be obtained as a 
limit of bounded solutions ($„, W„) of (|2.3ap . ()2.3bp in bounded domains Br^ with i?„ — )• oo as 
n — 7- oo. The reason for requiring the solutions to be obtained in such a way, is because we want 



these solutions to appear as a limit of bounded solutions of (|2.2ap . ()2.2bp as r — )• oo. Seemingly this 
implies that the mass at each aggregation point must be Svr. We would not attempt to give a precise 
meaning to these solutions in this paper, although it is likely that they could be given a precise 
meaning using some of the methods used in (H [T3l [TS] to define solutions of the two-dimensional 
Keller-Segel system for measures containing Dirac masses. Another alternative seems to be to use 
ideas analogous to the ones obtained in (7j. 

The solutions obtained in this paper will behave asymptotically as in ()2.4p . ()2.5p as r — t- oo. A 
particular case of these solutions corresponds to the case of radially symmetric solutions considered 
in [HI [9] . An alternative way of deriving the asymptotics of these solutions can be found in [19] . In 
this radially symmetric case, the corresponding solution of (j2.3ap . (|2.3bp has the form: 

<^r,s (y) = 87r6 (y) . (2.6) 

As was seen in [HJ [19], the solutions of (|l.lap - ()l.lbp with the asymptotics near the blow-up 
characterized by (j2.6p have the mass concentrated in a region of size: 



e(T) = /^e"V2, (2.7) 

2+7 

where K = 2e 2 with classical Euler's constant 7. The region where the mass aggregates can be 
described by means of a rescaling with a factor e (r) of the following stationary solution found in 

Us{x) = - -^, t;,(x) = -21og(l + |xn. (2.8) 

In this paper we will give most of the details concerning the asymptotics of solutions of (jl.lap - 
(jl.lbp which are bounded for t < T and blows up at t = T in the particular case of a limit function 
$s, a solution of (j2.3ap . ()2.3bp with the form (j2.4p concentrated in two peaks, (i.e. N = 2). The 
reason is twofold. First, the computations become more cumbersome for an increasing number 
of peaks, but without requiring essentially different ideas. On the other hand, the construction 
requires to check a sign condition for a suitable elliptic problem, as indicated in Remark II. 2 1 and 
this is what we have made numerically only in the case of two peaks. In any case, solutions of 
(|2.3ap . (|2.3bj) with the form (12. 4j) will be discussed in Section [H 

Due to the symmetry of the problem under rotations we can restrict ourselves to the case in 
which $s is given by: 

^s{y) = 87T[6{y-yi) + 6{y-y2)], yi = a, ya = -a, a =(2,0). (2.9) 

The detailed structure near the points y£, i = 1,2, can be computed by introducing boundary layers 
having many similarities to the ones described in [8] [19]. The rescaling factor e (r) will have a form 
similar to the one given in ()2.7p . although the value of the constant K will differ in general from 
the one obtained for the radially symmetric case. Actually, in the case of the asymptotics given 
by (|2.4p . the value of this constant could be different for each of the aggregation points. This will 
not be the case if $ {y, r) approaches the singular stationary solution $<; in (j2.9p due to symmetry 
considerations. 

A large portion of this paper consists in the detailed description of the boundary layers describing 
the regions of mass aggregation near the points 2/1,2/2- The computation of these layers will be made 
using the methods developed in [19] because the validity of some of the arguments in [8] is restricted 
to the radially symmetric case. 



We now describe shortly our strategy to compute the asymptotics of the solutions near the 
blow-up points. We will obtain outer and inner expansions for the solutions. The outer expansion 
is valid in the region where \y\ ~ 1 and \y — ye\ ^ e~"^ as r — t- oo, £ = 1,2, for some a > 
to be revealed later. The inner expansion is valid in the regions where \y — y(\ ~ e~"^, £ = 1,2. 
Both expansions are obtained under the assumption that the mass aggregating near the points yi 
concentrates in a region with width e^ (r) <^ 1, whose precise value will be computed later. Such 
assumption will be shown to be self-consistent with the derived asymptotics. There is a common 
region of validity where both outer and inner expansions make sense. The matching condition 
between both types of expansion in that intermediate region provides a set of differential equations 
for the functions £i (r) and these equations yield the asymptotics of such functions. 

We make extensive use of the asymptotic notation. We write / ^ (7 as x — )■ xq to indicate 
\\vD.x^xo f /g = 0, whereas / ~ (7 as x — )• xq to denote \\ra.x^xo f / 9 = 1- The notation f ~ g 
as X — )• xq indicates that the terms / and g have a comparable order of magnitude, that is, the 
existence of C > such that 1/C < lim infa;_5.2;Q f /g < lim sup^..^.^.^^ f /g < C. 

3 Inner expansions. 

3.1 Expansion of the solutions. 

We compute the asymptotics of the functions <I>, W defined in (12. lap . (I2.1b[) . In the case of radially 
symmetric solutions it is assumed that V<^ [yi, r) = with yi = 0. However, due to the lack of 
symmetry, points where the maximum of ^ are attained could change in time. We assume the 
existence of functions {y^ (t) : i = 1, 2, ..., A^} such that: 

V$fe(r),r) = 0, (3.1) 

lim ye (r) = y^. (3.2) 

It will be checked later that all these assumptions are self-consistent as usual in matched asymp- 
totics. Let us introduce the following set of variables to describe the inner solutions near each point 

ye ■■ 

C = ^^f^, (3.3a) 



1 

{edr)r 



^{y,r) = ^——^U{i,T). (3.3b) 



On the other hand, we will write, with a little abuse of notation, W {y^r) = VF(^, r). Notice 
that the variables .^, U (.^, r), and W {£,, r) depend on £, but these dependencies will not be explicitly 
written unless needed. Using ()2.2ap . (|2.2b|) . (I3.3a|) . and (|3.3b|) we obtain: 

ej^ = A^U - V5 (UV^W) + (2e,e,,, - ej) (u + ^^^ + (e,y,,, 
= A^W + U. 
We will now assume that the function ee (r) satisfies: 

Ei (r) ^ 1 as r — )■ 00, 

|q ^^I -|- \eiT-\ ■^ £e as r — 7- 00. 



f) ^^^' 


(3.4a) 




(3.4b) 




(3.5) 




(3.6) 



Assumptions similar to (j3.5p . (j3.6p are made in [19]. In addition, we will also assume in this paper: 

|y£,r| ^1 as T — )■ oo. (3.7) 

We now define in a precise manner the functions ee (r). We expect U, W to behave like the sta- 
tionary solution (j2.8p . The steady states of (jl.lap . (|l.lbp can be defined up to rescaling. Therefore 
the functions ei (r) could be computed up to a rescaling factor. The assumption C/(^,r) — )• Us(^) 
as r — )• oo would prescribe uniquely the leading order asymptotics of q (t) . Moreover, we can 
prescribe uniquely the function U, imposing the normalization: 

;7(0,t)=8 (3.8) 

or, in an equivalent manner: 

^iydr),r) = —^. (3.9) 

We then look for solutions of the system (13.4ap . (I3.4bp with the form of the following expansions: 

U (e, r) = Us (0 + f/i (C, r) + U2 (e, r) + C/3 (C, r) + C/4 (C, r) + ..., (3.10) 

W (e, r) = vs iO + Vf^i it t) + W2 (e, r) + ^^3 (C, r) + W^ (C, r) + ..., (3.11) 

where {ug, Vs) are the stationary solution as in ()2.8p . Notice that the function Vg is prescribed up 
to the addition of an arbitrary constant, but this can be ignored due to the form of the system 
(jl.lap - (|l.lbp . On the other hand, it will be assumed, as in [19], the terms Ui^Wi contain terms 
whose order of magnitude is e^ and that the terms U2 , W2 contain terms whose order of magnitude 
is (e^) or siyi^r up to logarithmic corrections like jloge^l ,r'^ or similar ones. Such logarithmic 
corrections will arise from terms like ei^r/^e or similar ones. The notation introduced in |19j and 
used also in this paper consists in writing all these terms as ef (w.l.a) (with logarithmic accuracy). 
We will include in Ui, Wi also the terms whose order of magnitude is ££ {w.l.a). Therefore: 

{Ui,Wi) ^ e£ (w.l.a) as T — )• 00. (3.12) 

On the other hand we will include in U2,W2 also the terms whose order of magnitude is e^y^^r 
{w.l.a) . Therefore: 

{U2,W2) ~ e? + £iye,T (w.l.a) as r — ^ 00 (3.13) 

In a similar manner, including in (Us, W3) terms of order e|, efye^r and Siyj^ (w.l.a) and including 
in (Ui,Wi) terms of order e^, e'jyi.T, ^"ivir (w.l.a) we obtain: 

(U3,W3) ^ ef + ejyi^r + eeylr (w-l-a) as r -^ 00, (3.14) 

(U4, VF4) ^ ej + e^ye^r + ^jyj.T (w.l.a) as r -)■ 00. (3.15) 

Making the assumptions (l3T^ - (l3TB it follows that the functions (Ui,Wi), (U2, W2), (U3, W3), 
and (C/4, VF4) satisfy respectively the following systems: 

= A5C/1 - V^ (ugV^Wi) - Vg (UiV^Vs) - ^V^us, (3.16a) 

= A^Wi + Ui, (3.16b) 



+ {2eiei^r - ef) i Us H ^ j + em^r'^^Us - -y^V^C/i , (3.17a) 

= A^W2 + U2, (3.17b) 

= A^f/3 - V^ (usV^Ws) - V^ {UiV^W2) - Vg {U2V^Wi) - V^ (C/sV^t;,) 

+ (2qq,, - ef) ("f/i + ^^^] " ^'^ + em,rV^Ui - ^V^U2 , (3.18a) 

= A5VI/3 + C/3, (3.18b) 

= A^Ui - V^ {usV^Wi) - Vg (f/iV^l^a) - V^ (t/aV^H^a) - V^ (t/gV^t^i) - V^ {U^VirVs) 

+ [letet^r - ej) (u2 + ^^j^) " 4^ + e^^irV^f/s - ^V^f/s, (3.19a) 

= A^W4 + U4. (3.19b) 

3.2 Computation of (f/i, VTi) , {U2, W2) . 

Due to (13. 8p we must solve (I3.17ap - (l3.19bp with conditions: 

C/i(0,r) = 0, C/2(0,r) = 0, [73(0, r) = 0, UiiO,T)=0. (3.20) 

We can easily obtain an exact solution of (j3.16ap - ()3.16bp : 

C/i(C,t)=0, W,i^,T) = -'-^^. (3.21) 

In order to compute {U2, W2) we notice that due to the linearity of (j3.17ap . (|3.17bp we can split 
its solution as: 

U2 = U2,l + t/2,2 + U2,3, W2 = W2,l + W2,2 + W2,3, 

where {U2j,W2j), j = 1,2,3, solve respectively: 

= A^f/s,! - V^ {usV^W2,i) - Vg {U2^V^Vs) + {2eiei^r - ej) (u^ + ^^^\ , (3.22a) 

= A^W2,i + U2,i, (3.22b) 

= Ag[/2,2 - V5 {UsV^W2,2) - V5 iU2,2V^Vs) + SiVi^r^^Us, (3.22c) 

= Ai:W2,2 + U2,2, (3.22d) 

= A^U2,3 - V5 {usV^W2,3) - V^ iU2,3^iVs) , (3.22e) 

= Ai:W2,3 + U2,3- (3.22f) 

We will check later that the term yg^j- is of order (e^) (w.l.a) . Therefore C/2,2 will be of order 

(q) (w.l.a) . Notice that this means that the terms Uk do not have a dependence (ei) (w.l.a) . 



On the other hand, at a first glance the system for (t/2,3) ^^2,3) could seem a bit odd for the 
absence of source terms. Actually (C/2,3, 1^2,3) will be chosen as a solution of the problem ()3.22ep . 
()3.22fp which are smooth for bounded values of |^|, but P1^2,3 becomes unbounded as |^| — )• 00. The 
contribution of (C/^2,3i 1^2,3) will be required to obtain a matching with some quadratic terms of the 
outer expansion having the angular dependencies proportional to {cos (20) , sin (20)} and giving 
corrections of order e'j {w.l.a). A detailed analysis of the matching conditions for the terms with 
this order of magnitude shows that, after a suitable rotation of the coordinate system, we may 
assume that the angular dependencies of the term (t/2,3, ^^2,3) are proportional to cos {29). We will 
then assume this angular dependence in the following. 

Due to (I3.20p we must have: 



f/2,fe(0,r)=0. 



1,2,3. 



(3.23) 



The solution of (|3.22ap . (I3.22bp satisfying the first condition in (I3.23P was obtained in [19] (where 
a slightly different notation was used). This solution has the form: 



U2,i (e, r) = Q2,i (r, r) , T^2,i (C, r) = ^2,1 (r, r) 



ICI 



where: 



51 (r, T)=r 



dV2A 



dr 

r dr \ dr 



+ Q 



2,1 



with: 



51 {r, r) = {2eeee,r - ef) 



{l + r^f Jo 



t2 



log (1 + t) 



l + t 



dt. 



According to the formulas (3.26) and (3.27) in [19j . we have the following asymptotics: 



Q2,i{r,r) 
dV2,x 



dr 



{r,T) 



{2eiei. 
{2eiee,- 



-?) 



^f) 



+ 



(log 



\2\-\ 



log (r^) 



+ 



as r 



(logr)^ 



00, 



as r 



00. 



The solution of the system (j3.22cp , ()3.22dp is given by the following simple formula: 

U2,2{^,r) = 0, W2,2iC,T)=eeytr-C 



(3.24) 

(3.25a) 
(3.25b) 

(3.25c) 

(3.26) 
(3.27) 

(3.28) 



We now consider the function (f/2,3j ^^2,3)- As explained before, this function, which is un- 
bounded at infinity, is just a homogeneous solution of the linearized problem. It will be needed 
due to the effect of the other singular points at the point under consideration. More precisely, the 
function W due to the points placed near y^ with k ^ i gives a contribution as |^| — )• 00 that 
will be matched with the term W2,3. The angular dependence of this term is cos {29) and its size 
ej {w.l.a). Therefore we look for a solution (f/2,3, ^^2,3) of ()3.22ep . p?22fP with the form: 



C/2,3(e, r) = Q2,3 {r, r) cos {29) , 1^2,3 (C, r) = ^2,3 {r, r) cos {29) 



(3.29) 



The system IK2M . ^^2f\i then reads: 

ldfdQ2,3\ 4 dUsdV2,3^^ „ dVsdQ2,3 ^ ,_„ . 

-TT r^^ - ^Q2,3 - -1 7^- + ^UsQ2,3 - -1 TT— = Oj (3.30a) 

r or \ or J r^ dr or dr or 

1 d / dVo3\ 4 

rYr V'^J " ~'^''' + ^^'^ = °- ^^-^"^^ 

The smoothness of (C/2,3, ^^2,3) at the origin (cf. also ()3.22ep . ()3.22fp ) imphes: 

Q2,3(0,T) = y2,3(0,T) = 0. (3.31) 

It was seen in [19] (cf. Theorem 13.21 below) that the space of solutions of (j3.30ap . ()3.30bp is a four 
dimensional linear space spanned by the set of functions {('i/'fc, Wfc) : k = 1, 2, 3, 4}. (We remark that 
the notation {ipk^Vk) was used in [19] instead, but we modify it here to avoid repetitions). The 
condition (|3.3ip implies that: 

(<52,3, 1^2,3) = Ki (^i,a;i) + Ks{'tps,U3) 

for some constants Ki, K3 G M. If K3 7^ 0, the growth of (^^3,^3) as |^| — )• 00 would imply 
that ($,VF) are very large for \y\ of order one, and this would contradict the hypothesis that <& 
approaches the steady state in (|2.4p as r — )• 00. Therefore K3 = and {Q2,3, ^2,3) is given by: 

Q2,3(r,r) = ^^^(rV3), V^.,3 (r, r) = ^ (r^ + 3) , (3.32) 

where i?2,3 = -62,3 (r) G M. Actually i?2,3 can be expected to be a function of r changing slowly with 
respect to this variable. By this we mean that -62,3 (r) does not have a factor like e~'^'^ with k ^ 0. 
The precise value of i?2,3 will be obtained later by matching the inner and the outer expansions. 
It will turn out to be of order (ei) (w.l.a) . Finally, notice that the formulas (I3.32P have been 
obtained for functions with angular dependence cos (26) , but similar formulas could be obtained if 
the angular dependence is replaced by sin (29) . The resulting coefficients i?2,3 will be denoted for 
functions with such an angular dependence as .62^3. 

In the following arguments, several more variables -64,2 ; -64,2; C3 (00) , .. will appear. They have 
some dependence on r, but we will not write this dependence explicitly unless needed. 

We remark also that the solutions of (|3.22ep . ()3.22fp cannot contain any radial contribution 
with angular dependence cos 9. Indeed, arguing as in the derivation of (j3.32p and using the fact 
that it is always possible to add a constant to V, it follows that such a contribution would yield 
an additional term in C/2,3 with the form Ki (r^ — l) (r^ + l) + K2r (r^ + l) cos 9. However, 
a Ki j^ or K2 7^ there would be a contradiction to (|3.ip . (|3.8p . Similar arguments exclude 
angular dependences cos (£9) with i > 2, since they would imply large values for <I>, W in the outer 
region where \y\ is of order one. 

3.3 Computation of {U3, W3) . 

Since Ui = and -V5 {U2V^Wi) - 2~''eeyeV^U2 = by ([321]), the system dSTSa]), (I3.18bp reads: 

= A5C/3 - V5 (usV^Ws) - V5 (UsV^Vs) , 
= A5VF3 + U3. 

9 



This system is similar to (j3.22ep . ()3.22fp . In order to obtain the matching of these terms with the 
corresponding ones in the outer region, we need an angular dependence proportional to cos (30). 
This dependence is the only one consistent with the rate of growth of these functions required to 
obtain the right matching with the outer part. We then write: 

Us (e, t) = Q3 (r, r) cos (39) , W^ {^, r) = V3 (r, r) cos (39) , (3.34) 

where {r,9) is as before. The function (Qs, V3) fulfills: 

~7r ^^r~ ~ ~lQ3 - -1 — TT- + ^UsQs - -1 — TT- = 0, (3.35a) 

r or \ or I r^ dr or ar or 



1 d / dVs 

r dr \ dr J r 



r-(^l^)-l2^3 + Q3 = 0, (3.35b) 



and the conditions implied by the regularity properties of U3, W3: 

Q3{0,t) = V3{0,t) = 0. (3.36) 

Using the solutions of these equations obtained in |19[ Theorems 4.1-4.3] we have: 

Q3{r,r) = ^^^{2r^+4), 1/3 (r, r) = -J^ (2^^ + 4) , (3.37) 

for some B3 G M. As in the case of i?2.3) ^3 could have some slow (meaning non-exponential in r) 
dependence on r. More precisely, it will behave like e^ (w.l.a) . We have just written terms with 
angular dependence cos {30) , but there are similar terms with dependence sin (39) characterized by 
means of a coefficient B3. 

3.4 Computation of (f/4, W4) . 
Using (IXT9iD . (|3l9b]) and ([MI]): 

=A^U4 - V^ (usV^Wi) - V^ ([/2V5VF2) - V^ {U3V^Wi) - V^ {U^V^Vs) + 

+ {2etei,r - e?) {1J2 + ^^^) " ^'^ + ^m,r^ i^^ - ^V^ZJs, (3.38a) 

=I\^W^ + [/4. (3.38b) 

Using (j3.2ip and (j3.28p we observe that: 

-V5(C/3V5W^i)-^V5t/3 = 0, 
-Vg (C/2V5W2,2) + em^^Vt.\J2 = 0. 
Then (fS^Hip . (|3.38bp yields: 

=l^^Ui - V5 {UsV^Wi) - Vg {U2V(W2,l) - V5 {U2V(W2,3) - 

- V^ {U^V^Vs) + {2ei£,^r - ej) (u2 + ^^\ - ej^, (3.39a) 

=A^W4 + C/4. (3.39b) 
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It is now convenient to split C/4, VF4 as: 

C/4 = ^74,1 + ^74,2, W4 = VF4,1 + W4,2, 

where: 

=A5C/4,i - V5 (n,VgW^4,i) - V5 ([/2,iVgTy2,i) - 

- V^ {Ui,iV^Vs) + {2eeee^r - ej) (u2,i + ^Mm 
=A^W4^i + C/4,1, 



,dU2,l 



(3.40a) 
(3.40b) 



= A^Ua,2 - V^ {usSI^Wa,2) - Vg (C/4,2V^t;,) + ^4,2 (C, r) , 
= A^VF4,2 + C/4,2 



(3.41a) 
(3.41b) 



with 



^4,2 (e,r) = -V5 {U2flVl-W2,l) - V5 ([/2V5VF2,3) + 



+ (2e^e^,^ - e|) ^72,3 + 



eVgC/2,3 



e^ 



= ?|i. (3.410 



The system (j3.40ap . (]3.40bp is the same as (3.16)-(3.18) in |19j and the solution can be obtained 
as indicated in that paper (although a slightly different notation for the functions is used there). 
The relevant information that we will need in this paper is the asymptotics of the solutions for 
large values of |^| which can be computed as follows. We define: 

l ^g2 5W4,1 _ 52 

r dr ^ dr r 

Then: 

„2 



u. 



4,1 



(3.42) 



92 — ££ (2efQ,r - ^£)t- 



r^logr 7r2 / „ 



2\2 



+ (2e£e^,T — ^e) 



\- O [ (log r 



(3.43) 
as r —7- 00. Similar asymptotic formulas can be obtained for dg2/dr. 

In order to solve ()3.41bp - (|3.41c|) we need to compute 5*4^2 Hi't)- Using (|3.24p and (|3.29p we 
obtain, after some elementary but tedious computations: 



where: 



^4,2 (e, r) = Gi (r, r) + G2 (r, r) cos [29) + G3 (r, r) cos (40) , 



Gi(rr)=--— frQ2 3^^ 
2r dr \ ' dr 

^ t . 4Q2,i^2,3 ^d ( n ^^^A ^ ^ f n ^^2,1 
r^ r or \ or J r or \ or 

,2^^2,3,.^ ^ ,2N 48i?2,3r^ 

dV2fl 



J.2 2r dr \ ' dr 



(3.44) 
(3.45a) 

(3.45b) 
(3.45c) 
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The form of Si^2 {^,'t) in p.44p suggests to split (^74,2,^^4,2) as: 

U4,2 = f/4,2,1 + ^4,2,2 + ^4,2,3, ^^4,2 = W^4,2,l + W^4,2,2 + ^^4,2,3 

with {(C^4,2,fc5 W^4,2,fc) : A; = 1, 2, 3} having the angular dependencies cos (2 (A; — 1) 9). Then: 

= A5C/4,2,fc - V^ (n,VgT4^4,2,fc) - Vg (C/4,2,fcV5^;,) + Gt cos (2 {k -1)9), (3.46a) 

= A^Wi,2,k + Ui,2,k (3.46b) 

with boundary conditions: 

C^4,2,fc(0,T)=0, k = 1,2,3. (3.47) 

The boundary conditions for 11^^2,2, ^4,2,3 are just consequences of the angular dependence of these 
functions and their smoothness properties, whereas condition (|3.47p for C/4,2,1 is just a consequence 
of (j3.8p . On the other hand, the angular dependencies of the functions Wi^2,2, W^4,2,3 yield: 



W^4,2,2 (0, r) = 1^4,2,3 (0, r) = 0, (3.48) 

whereas ()3.ip implies: 

^^^^^''(0,r) = 0. (3.49) 



dr 



3.5 Computation of (f/4,2,1, M^4,2,i) 



Lemma 3.1 Under the conditions ()3.47p and ()3.49p the system (j3.46ap . ()3.46bp mt/i fc = 1 /las a 
unique exact solution: 

o r^ (r'' + 4r2 + 9) 9W4 2 1 , ^or^(r^ + S) 

U,» = 2(B,,f^^^-^, _i^ = -(B,3)^_L_/. (3.50) 

where r = |^| and i?2,3 is the parameter in ()3.32p . 
Proof. Using (|3.45ap and (|3.46ap we obtain: 

= -T- r — ^ — rUs 7^^^ TT '^^4,2,1^- - Tj-TT ^^2,3 



r 3r \ dr J r dr \ dr J r dr \ ' ' dr J 2r dr \ ' dr 

Integrating this equation and using (I3.32P as well as (I3.47P we obtain: 

"^ dr ""'V dr "^^^'^'^ dr - 2 dr ' ^'^'^^^ 



On the other hand, defining 
and using (I3.46bj) . we obtain: 



^4,2,1 = '^ — ^^-^ (3.52) 

f/4,2,r = -^^. (3.53) 

The smoothness of the function ^4,2,1 implies -M4,2,i (0, r) = 0. Using (j3.47p we then obtain: 

M4,2,i (r, t) = o (r^) as r ^ 0. (3.54) 
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Plugging (1332)1 and ([333]) into (l33T]l we obtain: 

or \r or J (1 + 7^2^^ ' ^ i j^ j-^ Q^ 2 or 

In order to solve this equation we use the change of variables: 

M4,2, = ^^. (3.56) 

Note that (f33ij) imphes 

F4,2,i = 0(1) as r ^0. (3.57) 

Plugging (j3.56p into (|3.55p . we obtain: 



r' 5^^4,2,1 r . 2 . 9F4,2,i _ rQ2,3 9^2 



o\ ^-^4,2,1 ^ ' ^2,3 "■'Z.j .„ ^^x 

(1 + ^2)2 9^2 (^2 + 1)3 V -^^ 9^ 2 9r • ^ ^ ^ 

Using now (I3.32P to compute the right-hand side of (I3.58P we arrive at: 

52^^421 r2- 3(9^4 21 , .2 r^ [r'^ + ^) , A ■? ^ 

"■^TT^ - 2 ^ 1 n = -8 ^2,3' , ; J r-4 + 2r2 + 3 . 
or^ r^ + 1 or (7-2 _^ i)'^ ^ ' 

This is a first order linear differential equation for dF4^2,i/dr that can be integrated explicitly. After 
some computations we obtain: 

m,2,i _ 4(i?2.3)'r3(r^ + 3r2 + 3) 

dr (r2 + l)2 ■ ^ ' ' 

In the derivation of p.59p we have used that, due to (j3.57p . the value of -^4,2,1 (0, r) must be finite. 
Integrating now (I3.59P and using also (I3.57P we obtain: 



-^4,2,1 — — (-82,3) 



r^ fr2 + 3) 



r2 + l 
Using now ()3.54p and p.56p we have: 



M4,2,l = - (^2,3)' 



^6 ^^2 _^ 3) 



whence p.52p and p.53p yield the desired result. ■ 

3.6 Computation of (f/4,2,25 ^^4,2,2) • 

3.6.1 Reduction of the problem to ODEs for (54,2,2, ^,2,2- 

The functions 1/4^2,2, ^4,2, 2 satisfy the system (|3.46ap . ()3.46bp with k = 2 together with conditions 
(|3.47p . p.48p . In order to remove angular dependence we look for solutions of these equations in 
the form: 

^74,2,2 = Q4,2,2 {r, t) cos {29) , H^4,2,2 = ^4,2,2 (r, t) COS (26*) , 

13 



where (r, 6) is as before. It then fohows from (j3.46ap and ()3.46bp with fc = 2 as well as (j2.8p that: 

= ~7r ^^ ~ ^<34,2,2 + —: -3 ^ ^ 2 I 1 ^ ^ 7; 7:2'34,2,2 + G2, 

r or \ or J r^ (r^ + 1) ar r^ + 1 ar (l -)- j-^) 

(3.60a) 

0= -^- h^ ^ ^V4,2,2 + Q4,2,2- (3.60b) 



The precise formula of G2 may be computed, using (I3.25bp . (I3.25ap and (|3.45bp . as: 



8^2,3 ^- (r^ + 4r2 + 9) dgx 32^2,3 (r^ - 3) 
(r2 + \f dr ^ (r2 + 1)^ 



G2 = — ^^^^tV-:t3 — -4t + ;':\.4 ' gi+ 



8^2 3 (2e£e£ r - ef ) r-% . . x . ^ 8e?r2 (r^ + 3) , 

+ ,2'';,4 (r^ + 2r^ + 9) - (i?2,3). r 2 ^ ^3 ' ^^'^^^ 

where gi is the function as in (|3.25cp . The derivation of (|3.6ip requires just a long but elementary 
computation. On the other hand, the conditions (j3.47p and (|3.48p imply: 

Q4,2,2(0,r) = 1/4,2,2 (0,r)=0. (3.62) 

The system (I3.60ap . (|3.60bp is a nonhomogeneous linear system with source term G2 as in (I3.6ip . 
In order to study the asymptotics of their solutions we examine in detail the solutions of the 
homogeneous part of this system. 

3.6.2 Study of the homogeneous system. 

The homogeneous part of the system (j3.60ap . ()3.60bp can be written as the linear system: 

1 d / d^\ L^ ^ , 32r duj ^ Ar d^P ^ 16 , 



2 



^=Z-\r—]--^^ + i^ (3.63b) 



1 d / duj\ L 
r dr \ dr J r 

with L = 2. This system was studied in detail in |19j . Four linearly independent solutions 
{ipk,^k) 1 k = 1,2,3,4, were obtained and their asymptotics for large and small r were computed 
there. We need to compute an error term in the asymptotics of UJ4 for L = 2, 3,4, ... in a manner 
more detailed than in |19] . The following result is basically a reformulation of |19| Theorem 4.3]. 

Theorem 3.2 Suppose that L = 2,3,4,... A general solution of (I3.63ap . (|3.63bp is a linear com- 
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bination of four particular functions {{'tpkj^k) •> ^ = 1? 2, 3, 4}, whose asymptotics are given by: 

^i(r)= [(L-l)rVL + l], a;i (r) = -^-- [(L - 1) r^ + L + l] , (3.64a) 
(r^ + Ij r + i 

^^2(0= ^^ „^ [(L + l)r^ + L-l], a;2(r) = ^-4— T[(L + l)r2 + L-l], (3.64b) 

V'3 (r) ~ 8r^ as r -^ 0+, tJs (r) ~ -r^ as r -^ O'^, (3.64c) 

V'g (r) ~ 16ErLr^''^+^-2 as r ^ oo, wg (r) AKLr'^^^^ as r ^ oo, (3.64d) 

V'4(r)~8r~^ asr-;>0+, ^^4 (r) ~ -r^^ as r -;> 0+, (3.64e) 
V'4 (r) ~ 16CLr-^'^+^"2 as r ^ 00, a;4 (r) ~ Cl (^Lr-^ - Ar-"^^^^] + o {r'^-^) 



as r —)• 00 



(3.64f) 

for some real numbers Cl,Ki and k^. 

Remark 3.3 We have Cl,K£ > 0, whereas kl could be zero for some L. 

Remark 3.4 The only difference between Theorem \3.S\ and |ig|. Theorem 4-3J is that the last 
formula in ()3.64fp is written as 004 = o (r~^) as r —t- 00 in \19[ Theorem 4-3] ■ There is a typo 
there. The correct formula intended in that paper is u^ = O {r~ ) as r — )■ 00. Formula (|3.64fp 
provides a precise asymptotics for 0)4^. 

Proof. We need to prove only ()3.64fp . To show it we define, given a solution (^,w) of (|3.63ap . 
()3.63bp , two functions F, G by means of: 

^ = ^-^(F + G), u; = -^(F-G). (3.65) 

Then: 

(fG_ ( 1-2L 8r \ dG ( SL + 12 16 ^^_^ 

~d^ + \—r T^Tl) 1^ ^ [y^TT ~ {r-^ + lf) ' ^ ^ ^ 

d'^F (1-2L Ar \ dF 4 (L + 1) ^ 4 f dG ,^ , \ ^, , , ^, 

It was proven in [19] that there exists a unique solution of (j3.66p satisfying: 

Gp (0) = 1, Gp (r) ~ Clv^^ asr^oo, /3l=4 + L- ^/A + L^. 

Moreover, since the point r = 00 is a regular singular point for ()3.66p . we can use Frobenius theory 
to compute a power series expansion for G^ (r) as r — )■ 00 to observe: 



G^ (r) = Gl 



rh + 2VZ^T4-1 ^ 



^/IJTl+l 



h-2 + Q 1^^/3,-4^ 



as r — )■ 00 (3.68) 



as well as similar bounds for the derivatives. Two independent solutions of the homogeneous 
equation associated to (|3.67p are given by: 

Fi,fc (r) = (L + 1) r^ + (L - 1) , F2,,, (r) = r^^ [(L -l)r^ + {L + 1)] . (3.69) 
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We then look for solutions of ()3.67|) with the form: 

Fp (r) = ai (r) F^^h (r) + 02 (r) F^^h (0 

under the constraints on ai and 02: 

«'i (r) Fi,;, (r) + a'2 (r) F2,/, (r) = 0, 
a\{r)F'{r) + a'^{r)F'{r)=S{r). 



(3.70) 



AfM 



We set 

Fi,h{r) F2,h{r) 

Using then Cramer's formula as well as (j3.69p and (j3.7ip . we obtain: 



2(L + l)L(L-l)(r2 + l) 



2^2L-l^ 



a'l (0 



F2,h {r) 
' A^ (r) 



S{r) 



'^^^^-W)'^'^ 



r[(L-l)r2 + (L + l)] 

2(L + l)L(L-l)(r2 + l)' 

(L + l)r2 + (L-l) 



>S{r), 



2{L + I) L {L - I) {r^ + ly r 



2^2L 



-5(r) 



(3.71) 

(3.72) 
(3.73) 



Notice that r = is also a regular singular point for (j3.66p . Then Gp (r) = 1 + O (r^) as r — )• as 
well as similar estimates for the derivatives. We then observe that the function 5 in (j3.67p satisfies: 

\S {r)\ < C {l + r^^-'^Y r>0. (3.74) 

Note that /3l S (2,4) for each L = 2,3,... We can then obtain a particular solution of (j3.67p . 



choosing oi, 02 in (13.70p as: 
ai (r) = 



1 



a2[r) 



[(L-l)r?2 + (L + l)]7? 
2{L + l)L{L-l)J^ (7?2 + l)2 

1 r^[{L + l)rf + {L-l)\r] 



2{L + l)L{L-l) 



(7?2 + l)2,^2L 



S {t]) di], 
S (7/) drj. 



The convergence of both integrals is a consequence of ()3.74p . Notice that (j3.68p as well as the 
definition of 5 in (13.67P implies the asymptotics: 



S{r) 



'2 + ^/4TX2 



+ O (r^^ ^ ) as r — ^ 00. 



(3.75) 



Then: 



[(L-l).^ + (L + l)].^^, ML - 1) L^Cr^^-^ ^ ,,„_,. ^_,^^^ 



(r2 + 1) 



whence: 



ai(r) 



Cl 



2 + V4TX2 
[2 + L) + ^/4TX2 



j.l^L-2 



2(L + l)(2 + ^/4TX2 



+ / W{r])dr], 
Jo 
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with 



\W{r)\ < 



C 



(1 + ^5-/3^) 

for some C > 0. Since 13^ G (2,4) it then follows that f^ \W (r)| dr < oo. Therefore: 



ai (r) 



Cl 



(2 + L) + ^/4TT2 



2 (L + 1) (^2 + \/4Tl2 
as r — 7- oo for some kl G M. We then have: 



^''~' C,., 



L+l 



W {rj) drj 



ai{r) 



Cl 



[2 + L) + ^/I+l? 






2 {L + l) [2 + ^/1+1? 
On the other hand, using (j3.75p we obtain: 



+ 



L + l 



+ 0(r 



as r — )• oo. 



(3.76) 



[(L + 1) r^ + (L - 1)] r ^ ^^^ 4L^ (L + 1) C.r^^-3-2L ^ ^ _^^ 



(r2 + l)2r2i 
as r — )• oo. Therefore: 



2 + \/4Tl2 



a2{r) 



Cl 



(2 - L) + ^/4TX2 



„/3i-2~2L 



+ 0(r 



,/3i-4-2L 



as r — )• oo. 



(3.77) 



2{L-1)(2 + VA + 1? 
Combining (I3.68p - (j3.70j) . (13.76J) with (I3.77P we arrive (after some a bit lengthy computation) at: 



Fp{r) = CL 



A + KLV'' - ^^^±^r/^^"2 ^ Q (1) 



\/Z2T4 + i 

Using (j3.68|) and (|3.78|) in (j3.65p we obtain ()3.64fp . This concludes the proof. 



as r — )■ oo. 



(3.78) 



Proposition 3.5 Let Cl and Kl with L >2 be the constants as in ()3.64dp and ()3.64fp respectively. 
Then the following identity holds: 

ClKl = ,., , (3.79) 



Proof. Set 



M{r) 



VZ^TI 



/ Vl (r) ^^2 (?^) V'S (J^) "04 (?') \ 

wi (r) tJ2 (r) wa (r) a;4 (r) 

^'i(o v^2(o V'3(o v^ao 

\ ^ (^) ^ ('^) W^ (^) ^4 ('^) / 



Ai (r) = det {M (r)) . 



Using the system ()3.63a|) . (I3.63bl) we obtain: 

d/^L 



dr 



(r) 



2 4r 

- + 



r r^ + 1 



AL(r) 
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The solution of this equation is given by: 

A^ (r) = ^^ 

for some El G R. On the other hand, using ()3.64ap . (]3.64bp . (j3.64cp . and ()3.64ep we obtain: 
, , , 2iO(L + l)(L-l)L2 

Therefore E^l = -2^0 (L + 1) (L - 1) L^ and 

2iO(L + l)(L-l)L2 
^^ = 2f2_^u^ ■ (3.80) 

On the other hand the asymptotics of ipi, Ui, i = 1,2,3,4, given in (|3.64ap . ()3.64bp . ()3.64dp . and 
^^I^ imply: 

Al (r) = -2i°L (L - 1) (L + 1) VL^ + ^KlClt"^ (1 + o (1)) as r ^ oo, 

whence: 

El = -2^° (L + 1) (L - 1) l2 = -210l (L - 1) (L + 1) Vl2 + AClKl 

and the result follows. ■ 

3.6.3 Asymptotics of ((54,2,2, ^4,2, 2) • 

Lemma 3.6 For any fixed r, the problem ()3.60ap - (|3.60bp -( |3.62p has a one- dimensional family of 
solutions, parameterized by -64^2 = -64^2 (''")• Its asymptotics as r — )• 00 is given by: 

Qa,2,2 (r, r) = 16ir2i?4,2r2^-2 - 2v^C2i^2i?2,3 UeeEi^r - e] + '"-^^4) + O (-^1^ {w.l.a) , 

(3.81a) 

V^4,2,2 (r, r) = -4i?4,3i^2r2^ + ^^^ (2e,e,,. - e} + ^^^l) -' logr + ... {w.l.a) 

(3.81b) 

as r —7- 00, where C2 and K2 are the constants as in Theorem, \3.S\ and i?2,3 = -62,3 (r) is the 
parameter in (]3.32p . 

Proof. Although the functions (54,2,2 and V4^2,2 depend on r through the dependence on r of 
the source term G2 = G2{r, r) we do not write them explicitly in the proof, because the proof relies 
purely on standard ODE arguments and the dependence on r does not play any role. 
We look for solutions with the form: 

4 4 

(54,2,2 [r) ='^bi (r) Tpi (r) , 14,2,2 {r) ='^h (r) oji (r) , 

where the functions (ipijUi) are as in Theorem 13.21 with L = 2. We impose the constraints: 

4 4 

Y^b',{r)^P,{r)=0, Y.^[{r)u,{r) = 0. 

i=l i=l 



18 



Then: 

4 4 

Ql,2,2 (r) =Y,b^ (r) ^[ (r) , ^(2,2 (0 = E ^* (^) ^«' (^) (3-82) 

and using (j3.60ap . ()3.60bp as well as the fact that the functions {ipi^oji) solve the homogeneous 
system (j3.63ap . (|3.63bp we obtain: 



E y, (r) 4 (r) + G2 (r) =0, ^ ^^ (^) '^^ (^) = 0' 



i=l 



i=l 



We can rewrite p.82p . (j3.83p in the vector form: 

dB{r) 



M(r) 



dr 



S{r) 



(3.83) 



(3.84) 



where 



B{r) 



/ h{r)\ / \ 

b2 (r) .._ 

bsir) ' ^^"^^ -G2{r) ' 

\b^{r)J \ / 

and Ai (r) is the matrix appeared in the proof of Lemma 13. 5i We denote as Ai^ {r;S{r)) the 
matrix obtained by replacing the /c— column of M (r) by S (r). Cramer's formula then yields: 



^U^)^ ^"'l';;f'-'» . *^^l,2.3,4. 



(3.85) 



where A2(r) is the Wronskian given in (j3.80p with L = 2. In order to avoid lengthy formulas we 
will use the following notation. We will denote as I?2,m) m = 1,2, 3, 4, the following determinant: 



D2,m (r) 



ipi (r) ipj (r) ipk (r) 
Ui (r) ujj (r) u)k (r) 
uj'i (r) uj'j (r) w^ (r) 



, i,j,k£ {1,2,3,4} \{m}, i < j < k. 



Then: 



b'm (r) 



{-ir~''G2{r)D2,m{r)r^r' + l) 
3-212 



-, 771 = 1,2,3,4. 



(3.86) 



(3.87) 



Our next goal is to compute the asymptotics of the functions 6^ (r) as r — t- 0. To this end we 
first compute the asymptotics of -D2,m (r) and G2 {r) . Using (|3.64ap . (I3.64bp . (j3.64cp . (j3.64ep in 
Theorem 13.21 with L = 2 we obtain: 



^2,1 (r) 



(l+0(r)). 



£'2,2 (r) = -3 • 2*^ • r (1 + O (r)) , 

I?2,3(r) = -^(l + 0(r)), 
D2,4 (r) = -3 ■ 2^ ■ r {1 + O (r)) 
as r — )• 0. On the other hand, (|3.25cp and (j3.6ip imply: 

G2 (r) = 3 • 2=^ (3 {2eeee,r - ej) ^2,3 - (^2,3), ej) r^l + O (r)) 



(3.88a) 
(3.88b) 

(3.88c) 
(3.^ 



(3.89) 
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as r 



0. Combining (IHIHTll - (13:89]) we obtain: 

b'l {r) = -^ (3 {2e,ei,r - ej) i?2,3 - (i?2,3).ef) (1 + (r)) : 

b'2 {r) = ^ (3 (2££e£,r - ej) 62,3 - (^2,3), e|) (1 + O (r)) , 

3r 

^3 (0 = -7^ (3 (2e,e,,. - ej) i?2,3 - (i?2,3), ej) (1 + (r)) 
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3r^ 



as r 



^4 (0 = ^ (3 (2e^eAr " q') ^2,3 - (i?2,3). e') (1 + O (r)) 
0. Integrating the equations (I3.90ap - (l3.90dp and using (13.621) we obtain: 



bi (r) = A 



^ {2eeee,r - ej) B2,3 - {B2,3)^el 
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r2(l + 0(r)) 



62 (r) 



3(2e^e^,^-e|)B2,3-(52,3),e 



^».6 



24 



r^il + 0{r)), 



h3 (r) = /33 



3(3(2£,£,,.-e|)i?2,3-(i?2,3).£|) ^2 (1^0 
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(l + 0(r)) 



64 (r) 



3 (2e^Q^T- — e^) i?2,3 — (-^2,3)^ e 



^».6 
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r'' (1 + O (r)) 



(3.90a) 
(3.90b) 
(3.90c) 
(3.90d) 

(3.91a) 
(3.91b) 
(3.91c) 
(3.91d) 



as r — )• for some /3i, /33 E M to be determined. Notice that we can then compute the functions &„ 
by means of: 



hm. {r) = f3m + 



(-1) 



m+l 



3-212 



G2{v)r]'{v +V D2,m{v)dv, m = 1,2,3,4; 



/32 = /34 = 0. 



(3.92a) 
(3.92b) 



We now proceed to compute the asymptotics of the terms bm (r) as r — t- 00. To this end we need 
to derive the asymptotics of the determinants -D2,m- We begin with I?2,i- Expanding with respect 
to the first row of the determinant in (|3.86p we obtain: 



-02,1 {r) = V'2 {r) 



UJ3 (r) W4 (r) 
^3 i'^) ^i (r) 



4^3 (r) 



UJ2 (r) uji (r) 
^'2 i'^) ^i i'^) 



+ ^4 (r) 



UJ2 {r) UJ3 (r) 
^'2 i'^) ^3 i'^) 



Using (I3.64bp . (I3.64dp . and (f3:64fl) we obtain: 

^2^"^ J,{r) J,{r) -""V 

Mr) :;t\ "?!:! =-3-2^K2C2.-(V2 + l)+o(.-) 



(3.93) 

(3.94) 
(3.95) 



as r 



as r 



00. Using the asymptotics ()3.64bp and ()3.64fp we obtain: 

UJ2 (r) oj'^ (r) - 4 (r) 6^4 (r) = 2^ • 3(^/2 - l)C2r-2v^-3 (^ ^ ^ ^^^^ (3 gg) 

00, whence, using (j3.93p - (|3.96p and taking into account that 2v2 < 4, we conchide: 

^2,1 (r) = -2^ • 3C2i<:2r~^ (1 + o (1)) as r ^ 00. (3.97a) 
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The asymptotics of -D2,m (^) ) m = 2, 3, 4, can be computed by a similar manner. The following 
asymptotics then follow: 

D2,2 (r) = 2^C2K2K2r^'^-^ (1 + o (1)) as r ^ oo, (3.97b) 

D2,3 (r) = -2^ • 3C2r-2v^-3 (1 + o (1)) as r ^ oo, (3.97c) 

D2,4 (r) = -2^ • 3is:2r^^~^ (1 + (1)) as r ^ oo. (3.97d) 

Our next goal is to compute the asymptotics of the functions bm for large r. To this end, we first 
compute the asymptotics of G2 (r) as r — ?• 00. Using (j3.6ip as well as (j3.25cp we obtain: 

G2{r)^-^^({2eiee,r-ej) + ^-^^ej\ as r ^ 00. (3.98) 

Using now (I3:97a - (l3.97dp as weh as I^M^, (13.64bp . (I3.64dp . dSHg, IK92^ . and (I3.92bp we 
obtain the following asymptotics: 

C2K2B2;i ( I 2\ , (-^2,3)^ 2\ 1 , ^q ^^^ 

^1 (^) '" 2 — \ (^^^^^'^ ~ ^^) + ^ ^ Q ) log r as r -> 00, (3.99) 

C2K2K2B2;3, //„ 2^ , (^2,3)^ 2\ 2^2+2 , /o inn\ 

^"^^ ^ 3 . 24 ( V2 + 1) V ^ ^^'^ -^i)+ Q^ ^ eAr'^^ as r ^ oo. (3.100) 

In the case of 63 (r) we have that J^ \G2 {y])\ yf {rf + l) |-D2,i (^)l dry < 00. Assuming that ^-^ = 
O (^ef) {w.l.a) , as it corresponds to this class of terms we would then obtain: 

fe3 (0 ^ i?4,2 = /33 + ^7^ y G2 (r/) r?2 (ry2 + 1) 2^2,3 (r/) ^r/, (3.101) 

where i?4,2 = O (ef) {w.l.a) uniformly for large r. Finally: 

&4(r)~-^j^|^((2.,s,,-e?) + ^e?)r2v/^+^ as r ^ 00. (3.102) 

We now compute the asymptotics of ^4,2,2 {r) , 14,2,2 {r). We use (I3.64al) . (|3.64bp . (|3.64dp . (l3^64f)) 
combined with ()3.99p - ()3.102p to obtain the asymptotics of V4,2,2 as in (|3.81bp and 

n (r-\ -^(\K R ^2^2-2 (^2-^2^2,3 //„ 2N , [^^^^2\ 

<^4,2,2 i?^) — iOA2-D4^2^ p (^E^e^.r " ^l) H 7^ ^i - 

V2 + 1 V ' ^2,3 / 

^^3(0 P^ .^^^2/ 2 , ,^2n__.^^^^ , ^/^ ^ 



3-2 
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y^" G2 (r?)r?2(r?2 + 1)2 1)2,3 (r/)dr? + 0(^-^) («;./•«) (3.103) 



as r — )• oo. It is important to remark that in the computation of the asymptotics of V4,2,2 there is a 
cancellation of the leading order of 62 {t) ^2 {f) + 64 (r) lo^ (r) . We now estimate the integral term 
on the right-hand side of (j3.103p . The leading order of the integral is then computed as: 

j G2 (ji) r,^ [r,^ + 1)' Z?2,3 {ri) dr? ~ 2^ • 3^2^2,3 [{2ete,,r - ej) + -^^'j -^J^-^ 
as r — )• 00. Combining this formula with (|3.103p as well as ()3.64dp we obtain (|3.81ap . ■ 

Remark 3.7 It will be seen later in Subsection \5.5\ that in the outer region the terms ^^2^4,2'^^ 
would give a contribution for <l> of order i?4,2e£ \y ~ Vel ^^ ^^e outer variables. In order 

for this term to be smaller than ef we would need -64,2 = 0{e^ ). This implies basically that 

-64,2 is very small in this region. 
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3.7 Computation of (f/4,2,3, 1^4,2,3) • 

We now compute the functions C/4,2,3) W^4,2,3 which satisfy (j3.46ap . (]3.46bp with k = 3 together 
with boundary condition (j3.47p . (j3.48p . We can ignore the presence of homogeneous terms of the 
equations, because all such terms can be included in the parameters -62,3, -62,3 (cf. (j3.32p ). We 
can assume also that the angular dependence of the functions C/4,2,3) W^4,2,3 is cos (40) : 



f^4,2,3 = (34,2,3 {r, t) cos (40) , ^4,2,3 = ^4,2,3 {r, t) COS (40) 

Plugging (j3.104p into the system (I3.46ap . (|3.46bp with fc = 3, we obtain: 



r dr \ dr 



16 dus dV4,2,3 , _ „ dQi^2,3dvs , „ 

-2*^4,2,3 1 ^ r ^ltsV4,2,3 JTI TT + *-^3; 



dr dr 



dr dr 



„ I d f 9^^\ 16 

= -^- r ' T ^4,2,3 + Q4,2,3- 

r or \ or J r^ 
The conditions p.47p and ()3.48p respectively imply 

Q4,2,3(0,t)=0, 1/4,2,3(0, r)=0. 



(3.104) 

(3.105a) 
(3.105b) 

(3.105c) 



Lemma 3.8 The problem (J3.105ap - ()3.105b]) - (j3.105cp has a two-dimensional family of solutions 
parametrized by ci(oo), 03(00). Moreover, its asymptotics as r — ?• 00 is given by 



^4,2,3 (r, r) = IGK^r^^-^c^ (00) + (lAci (00) + VbC^K^ (^2,3)') + o (1) , 



2ci (00) 



^4,2,3 {r, t) = -AK4CS (00) r^v^ + 3ci (00) r^ + 



C4K4 (^2,3)' \/5' 
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(3.106a) 

+ o{r'^) (3.106b) 



as r — )• 00, where C4 and K^ are the constants as in Theorem \3.B and i?2,3 
parameter in (j3.32p . 



B2,3{t) is the 



Proof. For the same reason as in the proof of Lemma 13.61 we avoid writing explicitly the de- 
pendencies on r in the proof unless needed. The solutions of the homogeneous equations associated 
to (j3.105ap . ()3.105bp have been described in Theorem 13. 2[ We then look for solutions of (j3.105ap . 
(j3.105bp in the form: 

4 4 

Q4,2,3 (r) =^Ci (r) tpi (r) , 14,2,3 (r) = ^Ci (r) ^i (r) 



i=l 



4 = 1 



under the constraint: 

4 4 

i=l i=l 

Arguing as in the proof of Lemma 13.61 we then obtain: 

(-l)'=G3(r)A,fc(r) 



where: 



DA,m (r) 



4 (r) 



Tpi (r) Vi {r) ipk (r) 
UJi (r) ujj (r) uJk {r) 
uj'i (r) uj'j (r) uj'^ (r) 



A4(r) 
, i,j,A:G{l,2,3,4}\{m}, i<j<k 
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and where A4 (r) is the Wronskian given in (|3.80p with L = 4. Hence: 



Cm (r) 



(_l)™+V2(^2 + l)-^3(^)^^^^(^ 



214 . 3 . 5 

By Theorem 13.21 we obtain the fohowing asymptotics: 

2'^-3 



-, m= 1,2,3,4. 



-04,1 (r) ~ 
-^4,3 (r) ~ 



2^ •3-5 



1)4,2 (r) ~ -2^ • 5 • r^, 
1)4,4 (r) 2^^ • 3 • 5 • r^ 



as r — ;■ 0. Using now (I3.37P and (I3.45cp we have: 

G3 (r) ~ 2^ • 3 • (-62,3)^ r^ as r ^ 0. 
Combining (I3.107P with (I3.108ap . (I3.108bp we deduce the asymptotics: 

c^{r)r , C2 (r) ~ 2 (^2,3) r, 

5 

4 (r) ~ -6 (^2,3)' r, c'4 (r) ~ 6 (.62,3)' r^ 



(3.107) 

(3.108a) 
(3.108b) 

(3.109) 



as r — 7- 0. Using then the conditions (j3.105cp as weh as ()3.64bp and ()3.64fp . we arrive at: 

ci (r) = 71 + — 

C2 (?^) = 



214 . 3 • 5 JO 
1 



214.3-5 JO 
1 



C3 (r) = 73 + ^14 

C4 (r) = 



214 . 3 • 5 JO 
1 



)14 



3-5 







e'(e' + i) G3(0A,i(e)de, 
eMe' + i)'G3(e)i)4,2(0rfe, 

e2(e' + i)'G3 (0^4,3(0^^, 
e' (e' + i)'G3 (6^,4(0 '^e- 



(3.110a) 
(3.110b) 
(3.110c) 
(3.110d) 



We then proceed to compute the asymptotics of Cm (?') as r — )• 00. To this end, we compute the 
asymptotics of G3 (r) and D^^m (r) as r — )• 00. Using (|3.37p and (j3.45cp we obtain: 



Gsir) 



32 (^2,3) 



1 + 0^ 



as r — ;• 00. 



In order to compute the asymptotics of D^^i (r) as r —t- 00 we write: 



1)4,1 (r) = il)2 (r) 



^3 (r) ^4 (r) 
^3 (r) ^4 (f) 



i^sir) 



L02 (r) ^4 (r) 
^2 (^) ^4 ('^) 



+ V'4 (r) 



W2 (?') 1^3 {r) 

^'2 (^) ^ (^) 



We now have, using Theorem 13.21 the following asymptotics: 

^C^KiKi (2 + Vt\ r^^~^ (1 + (1)) 
-40i^4 (1 + Vs) r^^"^ (1 + o (1)) 



^3 (r) u;4 (r) 

^ (r) w^ (r) 

^2 (?^) ^3 (j) 

^'2 {r) ^ (r) 



(3.111) 
(3.112) 

(3.113a) 
(3.113b) 
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as r — )• oo. In the computation of the second term on the right of (j3.112p we must take into account 
the cancellations of the leading order term. Theorem 13.21 yields: 



UJ2 (r) uji (r) 
^2 {r) w^ (r) 



5r-4(l + 0(r-2)) 



C4(K4r"^-4r-2^ + 



as r 



U2 (r) (1^4 {r) 
^2 ir) J^ (r) 



20r^5 (1 + (r"2^) C4 ('-4K4r-5 + S^Sr-^v^-i _^ q (-^-7 

00. Using that \/4 + L^ — L < 2 as well as the fact that L = 4 we obtain: 

40 (^ - 2) 04^"^^"^ + O (r^iO) (3.114) 

as r — )• 00. The use of (j3.112p - (|3.114p as well as the asymptotics of ^j in Theorem 13.21 yield 

D4 1 (r) = -5 • 2^ • \/5C4i^4r^^ (1 + o (1)) (3.115a) 

as r — 7- 00. On the other hand, a direct computation using Theorem 13.21 gives: 

D4 2 (r) = 3 • 2'^C4ir4K4r2^^^ - 3 • 2*^ • ^flC^K^r (1 + o (1)) , (3.115b) 

D4 3 (r) = -5 • 3 • 2^C4r-2v^^3 (^ ^ ^ (^^^ ^ (3.115c) 

1)4 4 (r) = -5 • 3 • ilK^r^^-^ - 5 • 3 • 2«K4r2^"^ (1 + o (1)) (3.115d) 

as r — )• 00. We can now compute the asymptotics of the functions Cm (r) as r — )• 00. Since ci(r) 
and C3(r) are convergent to finite numbers as r — )■ 00, we may write them respectively as: 



1 /'OO 

ci (r) = ci (00) - ^,4 3_^ I e' (e' + 1) G3 (6 i?4,i (e) dc, 



C3 (r) = C3 (00) 



with 



ci (00) = 71 + 
C3 (00) = 73 + 



:)14 



1 



3-5 



eMe2 + i)-G3(6A,3(e)dc 



214.3.5 /_ e'(e' + i) ^3(61)4,1(0^?, 







1 



)14 



3-5 



e'(e' + i) ^3 (6^,3 (Ode 



(3.116a) 
(3.116b) 

(3.116c) 
(3.116d) 



Since (|3.11ip and (j3.115ap imply 



it follows from (j3.116ap that 



ci (r) = ci (00) + 



y/bCiKi {B 



•2,3]^ 1 



2-3-4 



(1 + (1)) as r — )• 00, 



whence: 

ci{r)ijji{r) = 3ci(oo)r + 



V^C^Ki {B2,zf 



+ 2ci(oo) 



r + o(r ) as r — )• 00. 



(3.117) 
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The full formulas of C2(r),C4(r) given in ()3.110bp . ()3.110dp as well as Theorem 13.21 show that 

C4K4L>4,4(^)" 



C2 (r) UJ2 (r) + C4 (r) ^4 (r) 
1 1 



214. 3 r4 



Jo 



A,2 (0 + 



de- 



-wsv^of^) re{e + ifGs{()D,AOdc+ 



+ 



Cj 1 
214 . 3 . 5 ;^ 



4r-2^+^ + o (^ J_^ ^ ^' e' (^2 + 1) 2 G3 (0 1)4,4 (0 de 



A quick check using (|3.11ip and Theoreni l3. 21 shows that I2 grow at most with the rate of 0{r'^^^ 4) 
as r — )• 00. On the other hand, it is readily seen by (|3.111|) and ()3.115dp that 

as r — ;■ 00. To compute /i we note that: 

D,2{r) + ^^^^^^^ = -2,-2^-V^C,K,r{l + o{l)) 
5 

as r — )• 00 due to (I3.115bp and (|3.115dp . Using this as well as p.llip . we obtain: 



J. = ^g'^^t,'"^-''^.^(l + o(l 



)) 

as r —7- 00. Summarizing, we have obtained the asymptotics: 

CiKi{B2,zf VI -I 



C2 (r) a;2 (r) + C4 (r) a;4 (r) 



12 



\/5 + l 



+ o(r2) 



as r — >• 00. 
Since 



1 roo 

C3 (r) = C3 (00) - ^,,3^ I e {e + 1) Ga (0 1)4,3 (6 d^, 
it follows from ()3.64dp and (|3.11ip that 

C3 (r) = C3 (00) + f^i^r^'-^'"-' (1 + « (1)) 



as r — >• 00, whence: 



C3 (r) ^3 (r) = -4K4C3 (00) r^v^ - ^^^^^^ (1 + « (1)) 

2 (Vo - Ij 



(3.118) 



(3.119) 



as r — 7- 00. We can then compute the whole asymptotics of ^4^2, 3 {f) as in ()3.106bp . By (j3.117p - 
(|3TT9]1 we conclude (l3.1U6ap . ■ 
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Remark 3.9 To see the contribution due to C3 (r) ^3 (r) we need to study the asymptotics o/Q4^2,3- 
Using (I3.64dj) we obtain the matching condition: 

Q4,2,3 {r, t) ~ I6C3 (00) i^4r^^^~^ as r ^ 00. 

This contribution would give terms of order e^" = e^~ ^ e^ in the self-similar region 

where \y — ye\ is of order one. Then the contribution of this term to <I> would be much larger than 
one unless C3 (00) is small as r —)■ 00. It will be seen in Subsection \5.5\ that 03(00) = 0{e^ ^ ). 

4 Outer expansions. 

In the analysis of inner expansions we have derived the asymptotics of the solution for a general 
number of peaks, but we will compute outer expansions only for the case of two peaks for the reason 
stated in the previous sections. The notation of the singularities has been denoted by {y^}^^, which 
solves (j2.5p . in the analysis of inner expansions. In the particular case where N = 2 we will write 
yi = a and 1/2 = —^ with |a| = 2. We may assume, without loss of generality, that a = (2,0). 

In this section we derive outer expansions for the solution of (jl.lap . (]l.lb|) . i.e. for regions where 
\y\ is of order one. To this end we argue as in the derivation of (3.40)-(3.48) in [19]. We look for 
expansions with the form: 

^y, r) = ejn{y) + eiei,rZ{y) + • • • , (4.1a) 

W{y,T) = Wo{y,T)+ejWi{y,T) + --- . (4.1b) 

Then, to the leading order, we obtain: 

-AWo = 0, y^ye 

with matching condition: 

Vj, Wo (y) ~ - r^^ as y ^ ±a, 

where a = (2,0). Then: 

V,W„fe) = -(li^ + l<^). (4.2) 

Neglecting the terms of order ef we obtain: 

2 \\y-a\^ |y + a|2y 

Plugging the expansion (|4.1a|) into (j4.3p we obtain the equations: 

L{n) = 0, y/±a, (4.4) 

L{Z) = 2n, y/±a, (4.5) 

where 

L(Si, = A«-!L|5+fl(!^ + lfe±^Vv«-a (4.6) 

2 \\y- a\^ \y + a\^ J 
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Due to (j2.8p and (j3.10p we should impose the following matching condition for Jl : 



^(y) 



|y=Fa| 



as y 



±a. 



(4.7) 



On the other hand, in order to obtain matchings in the regions where |3;| is of order one, we 
must assume that Q (y) increases at most algebraically as \y\ — ?• oo. 

The function Q cannot be computed in this case by means of a closed formula as in the radial 
case considered in |19j . Nevertheless, it is possible to prove that the problem (|4.4p - (j4.6p - ()4.7p define 
uniquely a function with the properties required to describe the leading asymptotics of $ in the 
outer region. More precisely, the following result holds. 

Lemma 4.1 Assume that |a| = 2. Then for every Di,D2 G M there exists a unique solution of 
(J13D, (jMl) satisfying: 



^(y) 



D, 



as y — )• a, Q (y) 



D, 



\y-a\' |y + a| 

1^ iy)\ ^ |y|™ for \y\ > 5 for some m > 0. 



as y 



Moreover, the asymptotics of 0, near the singular points ita are given by: 



n (y) ~ Di 
n (y) ~ D2 



1 



|y-a| 
1 



|y + a| 

where A = A {Di, D2) is a constant and: 

\2" 



^1 (y - a) + $2 (y - a) + ^3 (y - a) + ^ 



+ ^1 (y + a) - ^2 (y + a) + *3 (y + a) + ^ 



(4.8) 
(4.9) 

(4.10a) 
(4.10b) 



^1 (Y) = — 
^ ^ 16 



2 ^(a-y)^ 



\Y\ 



\Y\ 



, ^2iY) 



96|y|^ 






, ^3{Y) 



1 (a-n' 
256 |y|4 



Remark 4.2 Concerning the constant A in the asymptotics (I4.10ap . (|4.10bp . we have computed 
its value using the PDE solver "PDE tool box" from the Matlab package. We have observed that 
the numerical value of A is between —1 and —0.9 for Di = D2 = 8. The crucial fact is that A < 0. 
This negativity is a sufficient condition to ensure that a certain differential equation satisfied by ee 
has solutions approaching zero as r —)• 00 . (See Subsection [5l 



Remark 4.3 A result analogous to Lemma \4-1\ could be shown using similar methods for singular- 
ities of 0, at arbitrary number of points, or more precisely for operators with the form: 



N 






with aj G M^, j = 1, ..., N, slj ^ a^ for j ^ k satisfying: 



4(y - aj) 



■ vn-n 



(4.11) 



afc 
2 



N 

E 



i=l, JT^k 



4(afc 
, |afc - 



|2 ■ 



(4.12) 
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If (|4.12p does not hold, the form of the asymptotics (j4.10ap , ()4.10bp would contain additional 
terms with the homogeneity of l/\y — aj|^. Actually the condition |a| = 2 in Lemma \4-1\ is just the 
condition ()4.12p in the case of two peaks, i.e. N = 2. The functions {'^i}^^ that would appear in 
the study of the general case (I4.11|) have similar homogeneity properties to the ones described in 
Lemma \4-i\ but slightly different functional forms. 



Remark 4.4 A characteristic feature of the expansions (j4.10ap . ()4.10bp is the absence of logarith- 
mic terms in ^3. Very likely this property holds in general under the assumption (I4.12p for every 
integer N > 2. We prove, however, such absence of logarithmic terms just in the case of two peaks. 



Proof. In what follows, we denote by Yi = y — a, ri = |Yi| , Y2 = y + a., r2 = \Y2\ for nota- 
tional simplicity. The key point is to define suitable sub- and supersolutions having the expected 
asymptotics near the singular points. To this end we define auxiliary functions Wj as: 



Wj (y) 



1 + M/i (Yj) + i-iy^' ^2 (Yj) + M/3 (Yj) + a;,- (Y,) 



ri{Yj), i = l,2, 



(4.13) 



where r] (^) is a C°° cutoff function satisfying rj {^) = 1 for |^| < 1, rj {£,) = for |^| > 2, < 77 < 1 
and where the functions Uj (Yj) will be defined later. We construct a supersolution ri"*" and a 
subsolution Q~ of the form: 



n+ (y) = D^Wi (y) + D2W2 {y) + K, 
n- (y) = DiWi (y) + D2W2 (y) - K 



(4.14a) 
(4.14b) 



with a constant J^ > to be selected later. Some explicit but rather tedious computations yield: 



L -4+^l(>S) 



{a-Y.f 2{a-Y,f 1 






''A,) 



0^1 iYi) + 



(Fi-y2) , 4(a-y, 



2"' " ' r^r^ 



+ 



rrij) 



A 1 



-+ 



+ 



64r| 



16 

r{j). 



lap (n ■Y2)-2{a- Yi) (a • Fa) + 4 (o • Y 



.2 (Yi ■ Y2) 



(4.15) 



for j = 1, 2, where r (j) = 3 — j for j = 1, 2. In the derivation of these formulas we have used: 



1,, „ / M 4y.r, 






-0) ^ v^ly v"l. 



as r. 



0, j = l,2. 



This formula holds due to the assumption that |a| =2. In all these computations we often use: 



^r(j) - 16 = r, 



4(-l)^(a.y,), j = l,2. 



It follows from (|4.15p that: 



^^.^...^^^ 



3r^ 



(a-y,)' 



+ 



as Tj — )• 0, J = 1,2. 
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We now use: 



A (^2) (Y) + -i^ (y . V) (*2) (Y) = J^-^^ 



3\Yf - {a-Yf 



as well as the fact that the terms in L that are not A and -ttj^j [Y ■ V) yield only lower order 
contributions. Therefore, after some computations, it follows that: 



1 



L -,+^,{Y,) + {-iy^'^2{Y, 



i+i 






Using then that: 



A1'3 + 



4(y-v)^3 3(a-y)2 (a-y)^ 



\Y\ 



16\Y\ 



16|y| 



as well as the fact that the remaining part of L gives only lower order contributions, we obtain: 

L(^,+^i{Y,) + {-iy-^'^2{Y,) + ^s{Y,)]=J2f^k^^^+9AY,) J = 1,2 (4.16) 

\ J J fc=o l^il 

with Qj G U^ (^2 ( (— 1)'' aj j and some suitable /3fc G M. Using a separation variables argument 
we can construct functions tjj, j = 1,2, with the form: 



-.(^.) = E-^-^ 



y 



i = i,2. 



The constant k^ is selected in order that functions Uj may satisfy: 



y 



y 



fc+i 



(4.17) 



We then define the functions PF,-, j = 1,2, as in (jiT3|) . It follows from (jiT6]) and (jiTTl) that: 

L{W,) = fj{y), j = l,2 



with 



Ij\\L° 



< 00. Using the boundedness of fj as well as the fact that the functions W^ 



J' 



j = 1,2, are compactly supported, we observe that 0"*", n~ in ()4.14ap . ()4.14bp are respectively 
super- and subsolutions of ()4.4p in M^ \ {—a, a} if K is chosen sufficiently large. 
We now define a family of domains Dg fi as: 



Ds,B^ = BR{0)\[Bs{-a)UBs{sL)], < 6 < I, R>8. 
Let us consider the following family of boundary value problems: 
L {ns,R) =0 in Ds,R, 



(4.18) 



(4.19a) 



^S,R = Dj 



1 



+ ^1 (y) + (-1)^+1*2 {Yj) + ^3 (Yj) 



ns,R = on OBr (0) . 



on dBs [i-iy+U) , 3 = 1,2, (4.19b) 

(4.19c) 
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Classical results on elliptic equations (of. [6l Corollary 9.18] for instance) show that the func- 
tions ^s^Ft are uniquely defined for any 6 and R in ()4.18p . Moreover, since fi" < 0,s^r < ^2"*" on 

[Jj=i 2 ^^s { (—1)''''' a j U OBr (0) ior K > sufficiently large independent of 5 and R, it follows 
by comparison that: 

n- < ns,R < n+ m Ds,r. (4.20) 

Classical regularity theory for elliptic equations implies that |V'^il5^/j|, k = 1,2,3, are bounded in 
compact sets of Ds^r. A compactness argument then shows that there exists a smooth function il 
satisfying Q~ < il < Q,^ in M^ \ {—a, a} and a subsequence {{6e,R£)} such that (6i,Re) — )• (0, oo) 
and: 

^5e,Re -^ ^ as ^ — )> OO in C^ (/C) 

for each compact JC CM?\ {—a, a} . Therefore L (il) = in M^ \ {—a, a} . Moreover, the functions: 



cPj [y) =n{y)- D, 



^ + *i (Yj) + {-iy+'^2 {Y,) + ^3 {Y,] 



i = l,2. 



are bounded in a neighborhood of the points {—a, a} and satisfy: 
where 



Qj {y) 



\Y^ 



\Qj (y)| <C in < 1^^! < 1, j = 1,2. (4.21) 

To conclude the proof it only remains to show that the limits Imiy^a 4'1 (y) and liuiy^-a <i>2 (y) 
exist. To this end we estimate 
^R,j (.0 = (Pj (±a + RO ■ Then: 



exist. To this end we estimate the derivatives of <j)j as follows. For each < i? < 1 we define 



^i^R,j + ^^ " j^f""'' + aR (e) • y^^R,, = 0{R), 1<\C\< 4, 

where \aR (^)| < CR. Classical regularity theory for elliptic equations yields |V^9?i?| < C in 1/2 < 
1^1 < 2. Therefore: 

\(l),{y)\ + \Yj\\V(l){y)\<C in < |y,| < 1, i = l,2 (4.22) 

for some C > 0. In order to prove the existence of the limits limy_>.-|-a (pj iy) we now use a Fourier 
analysis argument. We use polar coordinates defined as: 

Yj = y ^ a = pj (cos 6j,sm 9j) , i = 1, 2. 

We then write: 

oo 
n=— oo 

The functions c„ (pj) solve a second order ODE, which can be solved explicitly: 

c..(ft)=^,,„p-*+^.„pf - n^f'^^d.- r (ft r-" 4-4..,, =1,2, 



Pj " ■' ' Q^n Oin J Pj ^ S / OLji C^r, 
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where: 

f-27r 



Q • „ (s) = J- f Q. (±a + s (cos 9i, sin 0,)) e^*"^^(i%, (4.23) 

27r Jo 

a+ = -2 + \/4 + n2, a" = -2 - y^IT^, n G Z (4.24) 

and where Ai^n, ^2,n are constants related to the Fourier coefficients of the functions ^j {l,9j). 
Since these functions are in C°° {'^^)j fo^ every /3 > 0, there exists a constant C/j > such that 

\Ai,n\ + \A2.n\ < — ^^, ri€Z. (4.25) 

l + |n|P 

On the other hand, due to (j4.2ip and (j4.22p the coefficients c„ (pj) and Qj^n {Pj) are bounded for 
< /9j < 1. This imphes: 



'0 V^y On - On 

whence: 

c„(«) = A,„P?--/'(^)"'%^<i.+ r(^)°'^^<!., ; = i,2, ue 

Jpj ^ s J aZ - an Jo ^ * ^ On - an 



Using (jOTIl - ifOS]) we obtain: 



,^5-2 „ n^V5-2 



, _ , C^Pr~ Cpi Cp^ 

\Cn (pj) - Ai^oKol < ^ I |2 +- 7^ < 



111 1^ 1 ll 1-^ 1 ll 1^ 

+ |n| 1 + |n| 1 + \n\ 

where symbol 6nfi stands for the Kronecker delta. Then: 

oo 

\ct> (y) - ^i,o5n,o| < Cpf-' Yl -r-Tl2 < Cpf-' 

and the therefore the limits limy^-i-a 4'j iv) exist. The fact that the value of A is the same in (|4.10ap 
and ()4.10bp follows by a symmetry argument. 

To prove the uniqueness result we construct a supersolution for (|4.4p . Consider a function $7"^ in 
(j4.14ap with K sufficiently large. We then modify the function $7+ so that the constant K becomes 
the polynomial \y\^ for large values of \y\. Since the main terms in the operator L for large values 
of \y\ are 2~^y ■ Vfi and — Q, the modified function Cf^ satisfies L(Q"'") < 0. This modification is 
possible, because these leading terms yield positive contributions. The difference of two solutions 
of (j4.4p satisfying (j4.8p may be estimated by eVf^ for y — )■ ita and for \y\ = R with e > arbitrarily 
small and i? > large enough. A comparison argument then shows that the difference is bounded 
by £0+ in the regions Br (0) \ B^ (ia) for 6 small. Taking the limit e — t- we know that both 
functions are the same, whence the uniqueness follows. ■ 

Remark 4.5 Equation (14. 4p suggests that Q{y) r^ ip (6) / \y\ as \y\ — )• oo, for some function ^ (6) 
whose precise formula does not seem easy to derive. However, we will not attempt to compute this 
function in detail in this paper. 
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We also need to study the function Z, a solution of (|4.5p . satisfying: 

1 



Z{y) = o[ -J as y ^ ±a, 

|-Z^(?/)| ^ lyl"^ for \y\ > 5 for some m > 0. 



(4.26) 

(4.27) 



Lemma 4.6 Suppose that |a| = 2. Let Q., Di, and D2 be as in Lemma \4-l\ Then there exists 
a unique solution of (|4.5p satisfying (|4.26p and (|4.27p . /is asymptotic behavior near the singular 
points {—a, a} is given by: 



Z{y)^Di 

Z {y) ~ D2 
for some constant B G 



[ 2|y-aP + 8^°s"^ "" 


1 (a.(y-a))2 " 
16 |,-a|^ +""] 


[ 2|, + aP + 8^°gl^ + ^l 


1 (a.(y + a))2 , „' 

i« , ,2 +-^ 
Id y + a 



as y 



as y 



(4.28a) 
(4.28b) 



Proof. The proof is similar to the one of Lemma [4. 11 Due to the linearity and the symmetry of 
the problem it is enough to consider the case Di = 1, D2 = 0. We can obtain sub- and supersolutions 
of (j4.5p with the help of the auxiliary function: 



W{y) 



11 1, „,, 1 1 {si-YiY 



2|>1| 



8 16 \Y^Y 



ri{Yi). 



(4.29) 



where rj (^) is a C°° cutoff function as in Lemma |4. II We then construct sub and supersolutions in 
the form: 

Z+{y) = W{y) + K, Z' (y) = W (y) - K. 

The terms between the brackets in (j4.29p have been chosen in order to balance terms of the function 
Q (y) . This requires some tedious, but otherwise straightforward computations. Arguing then as 
in the proof of Lemma |4. II we obtain the result. ■ 



We also need to study the asymptotics of the function Wi in (14.1bp . which solves the equation: 

- AWi = n, 2/ / ±a. (4.30) 

We obtain the following result: 

Lemma 4.7 Suppose that |a| = 2. Let Q, be as in Lemma \4- 1\ Then for every M^^^ ,M^^ G M, 
there exists at least one solution of (I4.30p satisfying: 

Wi (y) - DiG (y - a) - I^iMJ"^^ log \y-s.\=0 (1) as y ^ a, 



Wi 
(-a). 



Wi (y) - DiG (y + a) - DiM\ ^> log |y + a| = O (1 
hm MM = o, 



as y 



(4.31a) 
(4.31b) 

(4.31c) 
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where 

GC>') = --h?i - J (i°e I^D' + ^"s (29) (i-3id) 



,|^.4<'°''l^l>'+3i 



and where = 6 {Y) is the angle between the yi axis and Y. 

Moreover, two arbitrary solutions of (|4.30p satisfying (|4.31ap - (j4.31cp differ by a constant. We 
have the following asyniptotics for Wi (y) as y — t- a and y — t- —a : 

Wi (y) = D^G (y - a) + DiMI%^ log |y - a| - ^ \y - a| cos (30(a)) + 

+ aS") + Z^i/fj") . (y _ a) + O (|y - a|2) , (4.32a) 

Wi (y) = D2G (y + a) + i^zA/J^^J log |y + a| - ^ |y + a| cos (30(_,)) + 

+ a["'^ + IJsi^i"''^ • (y + a) + O (|y + aj^) , (4.32b) 



where Ofa)^ ^(-a) '^'^^ ^^^ angles between the horizontal axis and the vectors y— a andy+a respectively. 

r(a) ,f(-'^) 



(a) (—a) 9 (a) (—a) 

The vectors Kl , K^ G M and the constants A\ , A\ G M depend on an affine manner on 



the values of M^ -^ , M^ y^' . 

Proof. The proof is similar to the one of Lemma 14. li It reduces just to compute explicitly 
the solutions of the Poisson equation having as sources the terms in the asymptotics (I4.10ap . 
(|4.10bp . After removing the effect of these singular contributions, it only remains to obtain a 
solution of Poisson equation with a source term bounded as C/(l + \y\ ). This can be made 
using a supersolution behaving as C(log(|y|)) for large values of |y| . The uniqueness result is a 
consequence of Liouville's theorem for the Laplace equation. ■ 

5 Matching of the different terms 

In this Section we match the different terms in the inner and outer expansions and consequently 
derive evolution equations for the functions ££ (r) providing the width of the peaks. We will assume 
in the following that, due to symmetry considerations, all the functions £1 at the different peaks 
are the same. In general this does not need to be so. Moreover, there are non-symmetric singular 
self-similar solutions (cf. Section [2]) for which the corresponding values of the functions e^ (r) cannot 
be expected to be the same. The question of determining the relative sizes of the functions ££ (r) 
is interesting, but it will not be considered in this paper. Due to ()3.3bp this question is equivalent 
to determining the relative sizes of the maximum value of the function u at each of the different 
peaks (Notice however that all of them have the same mass Svr). 

We now describe how to match the different terms in the asymptotics as |.^| — )■ 00 of the 
expansions (|3.10p . (|3.1ip . We begin with the leading order terms. Since we restrict our analysis to 
the case of two peaks we assume in the following that ei = £2 = £■ We write for further reference 
the expansion of VyWo near y = a (cf. ()4.2p ): 

4y 2a Y 2a(a-y) jypa 2(a-y)^a {a-Y)Y 

\/ym - --rrp72 -7-^2-7-^2+ T^ + T7^^ i-[6 + — iZii 

II 1^1 1=^1 1^1 ^ 1^1 1^1 1^1 

(a-y)|y|2a |y|4a 2(a-y)^a |y|^y (a-y)^y 

|a| 1d|c^| 1=^1 4 |a| |a| 

where y = y — a and we have kept in this formula all the terms until third order in \Y\. 
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5.1 Leading terms. 

The leading order in (j3.10p . ()3.1ip is respectively given by the functions Us {(,) , Vg (S,)- We will 
denote as ^o,match, Wo,match the terms to be matched in the intermediate region |^| » l,|y — y^l ^ 
1 due to these terms in the expansion. Keeping just terms of order e| {w.l.a) in the region where 
\y — ye.\ becomes of order one, we then obtain: 

^0,matc/i (y,T) = ^, V yW^^rnatch (?/, t) = ■ — ^. (5.2) 

\y - yt\ \y - yA 

The matching of the term '^ yW^^match has been already taken into account in the derivation of ()4.2p 
that gives the asymptotics of the chemical field for |y| of order one up to corrections of order e|. 
On the other hand, due to (j4.10ap . ()4.10bp we obtain the matching of (|5.2p with (j4.1ap . assuming 
L»i = 8. 

5.2 Terms coming from Ui, Wi. 

We now match the terms C/i, Wi in ()3.10p . (j3.1ip with suitable terms in (j4.1ap . ()4.1bp . respec- 
tively. We denote as ^i^match-, ^i,match the terms to be matched in the intermediate region 
1^1 ^1, \y — ye\ ^ 1 due to these terms in the expansion. Notice that (j3.2ip shows: 

^l,match [y, t) = 0, VyWi^rnatch {y,V = -^■ 

We only need to match the term ^yWi^match with some of the terms in (|4.2p . Let lim-j-^ooyi = 
ye = a, since the case lim^_!.oo yi = —^ can be treated in a symmetric way. The most singular term 
of ()4.2p has been matched with VyWo^match- The next order in the expansion of V^Wo is —2a/ |a| 
(cf. (j5.ip ) and this matches with —yi/2 = —a/2 if we impose |a| = 2. Therefore matching of the 
terms of order e (w.l.a) in the region where |^| is of order one becomes possible if we impose that 
the drift terms due to the change to the self-similar variables and the chemotactic terms balance 
with each other. 

5.3 Terms coming from U2, W2. 

Let us denote as ^2,match, W2,match the terms appearing in the matching condition arising from 
the terms U2, W2 in the inner expansion. Using (13.26p . (j3.27p . (j3.28p . and (I3.32P we obtain the 
following formulas in the intermediate region ei <^ \y ~ IJil ^ 1 '■ 



^2,match (y, t) ~ {^EiEi^r " £?) 



^ +0' ' 



I - |2 \ I - |4 

\y-yi\ \\y-yA 



-2 



+ (5.3a) 



, 8^2,3 cos (20) 8^2,3 sin (20) 

\y - yi\ \y - yt\ 

W2,match [y, t) ~ ye^r ' [y " m) + V'2,1 H 5- \y - yi\'^ cos {26) -\ j- \y - yi\'^ sin {26) + ..., (5.3b) 

where V2,i is a radial term. It is more convenient to rewrite ()5.3bp in cartesian coordinates: 



B. 



2,3 



W^2,match (y, t) ~ yg^r ' (?/ " yi) + ^2,1 ^ f- (?/! " 2/£,l) " (2/2 " ^^2)' 

^ 1} 



\2 / - \2 



+ 



H ^{yi-yi,i){y2-yii,2) + - (5.4) 
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The last two terms of this formula can be matched with the quadratic terms of the expansion of 
VyWo(y) near the points y^. Using (15. ip it follows that VW2.match matches with Vj,Wo(y) if 



£2 



52,3 = -|, ^2,3 = 0, (5.5) 

where we use that |a| = 2. Using (I5.5P in (j5.3ap and transforming the resulting formula to cartesian 
coordinates we obtain: 



^2,match {y, t) ~ £(^ 



2,{yi-2f + {y2f 

\y - vA 



-i^ + of^^)... K/.a) (5.6) 

\y-yA \\y-yt 



and the first term in ()5.6p matches exactly with the term in the outer region multiplying \^i (y — a) 
in (|4.10ap due to the fact that Di = 8. 

It is illuminating to compute yi^-y in the first term of (j5.4p . matching the first term on the right- 
hand side of this formula with one of the terms in the outer expansion ()4.1bp . We will examine 
the case in which limT-_j>oo yt = a, since the case in which limT-_s.oo yt = —^ is similar. Using (I4.1bp . 
()4.32ap as well as the fact that Di = 8 we obtain the following terms in the outer expansion of 
'^yW which require to be matched with terms from the inner expansion: 



e} 



AY Y I y-L ^m\% Y 

^ - 2 (l°g 1^1) ^ - 2 '"^ (2^('^)) ^ + ^^+ 

n") Y sin ('36'/„^'l y-L 

(5.7) 



(^■. cosiW/^^) Y sin (36'(„i) y-L 
+ 8K;^) - -^-^ + -1-^ +0{\Y\) 



where Y = y — a. and we define Y = {—y2,yi) for Y = {yi,y2) ^M?. 

The first term in (j5.7p matches with a similar term coming from the function Vs in (|2.8p using 
Taylor series as |^| — t- oo. The second term in (j5.7p matches, to the leading order, with the first 
term on the right-hand side of (j3.27p . The third term in (j5.7p matches with the first corrective 
term that results in the Taylor expansion of l^^s in (|3.32p as |^| — )• oo. Notice that we use also in 
this matching (j5.5p . The matching of the term 8M| yy Y/ \Y\ plays a relevant role in determining 
yi^T-. Indeed, the contributions of similar order in the inner region are due to the terms log (r^) /r 
and — 2/r in (j3.27p . Due to the change of variables r = |,^| = \Y\ jei it follows that, to the leading 
order SMJyy = loge^. Lemma Wn\ thus yields K^ = of' loge^ as r — )• oo to the leading order. 
We can then match the term %K^ e^ in (|5.7p with the first term in (|5.4p . whence: 

/oo 
e\{s)\o^[Ei{s))ds asT->oo. 



This gives the desired asymptotic formula of the peaks stated in (j3.2p . The terms with the angular 
dependence W in ()5.7p are matched with some of the high order corrections coming from (j3.37p . 
However, this terms give smaller contributions and we do not pursue this computation in detail. 

5.4 Terms coming from [/a, VFs. 

We now match the terms coming from C/3, W3 which can be computed by means of (|3.34p . (j3.37p . 
We notice that, to the leading order, C/3 must match with the term e\^2 (^) in (|4.10ap . ()4.10bp . 
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Using that -Di^2 (^) = —6 ^ cos (36) / \Y\ (since Di = 8) we can match this term with the leading 
matching term coming from C/3, which can be written as (cf. (13.370 ): 

££ \y - a| 
whence: 

B3 = -^, Bs = 0. (5.8) 

We see that this gives also a matching for the terms in (|5.ip with angular dependence cos(30(a))! 
sin(30(a))- Using that |a| = 2 we can write those terms as: 

lypa 2(a.-Y?a. (a-Y)Y IFI^ 

' ' ^ ^ I V J _\ \ [(_ cos (20), sin (20))]. (5.9) 



2|a|" |a|6 |a|^ 2^ 

On the other hand, we can compute two terms of the asymptotics of Vg(V3(r) cos(30) as |^| — )■ 00 
using Taylor series. Rewriting the resulting expansion using the y— variable we obtain the following 
terms to be matched from the inner expansion: 

^^^r2(cos(20),-sin(20)) + — ( 2cos (30) ^ - 6sin (30) ^ ) . (5.10) 



ee ' ' '' ' " ee \ ' ' \Y\ ' ' \Y 

Using (|5.8p we obtain that the first term in (|5.10p matches with the term in (|5.9p and the second 
one matches with the terms in (|5.7p with angular dependence 30. 



5.5 Terms coming from ^4,^4 

We now match the asymptotics as |^| — t- 00 in the terms U{^,t), V {^,t) with the terms in 
the outer expansions (j4.1ap . ()4.1bp that are of order e| (w.l.a) as y — t- ita. These are terms in 
the outer expansion multiplying ^3 (Y) and A in (j4.10ap . ()4.10bp as well as the terms multiplying 
16^^ (a • (y — a)) / \y — a\ and B in (j4.28ap . (j4.28bp . Therefore, using also that Di = 8, we obtain 
that the outer expansion for <1> to be matched as y — t- a is: 

1^^^ + 8Aej - ''I''- (^ -J/ + 8Be,e,,r + e.e,,^ log |y|, (5.11) 

•^ I y I ^ \Y\ 

where y = y — a. 

Concerning the inner expansion we notice that the only radial terms giving contributions of 
order ej {w.l.a) in the matching region are the terms C/4^1 + Ki,2,i- Using (j3.42p . (j3.43p . and (j3.50p 
as well as the change of variables we obtain the following radial terms for $ to be matched: 

_2N /log|y| loge^ 5\ (2e(,e(,^r - ej) ej 



- (2.,.,, _ 4), ^^^^ - ^ - ^ j + A^r^^ig^ + ^, (5.12) 

where we have used (15. 5p . On the other hand, we can decompose the terms in (15. lip in radial terms 
and in terms with angular dependences cos (20) and cos (40). Using also that |a| = 2 we observe 
that the radial terms are: 

— f + ^Ae} + (85 - 1) EiEi^r + e^e£,r log \Y\. 
lb 
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We notice that the term containing log \Y\ can be matched, to the leading order, with a similar 
term in (I5.12p . On the other hand, the matching of the remaining terms provides an equation for 
e£ in the same manner as in I19j : 



^ '- ^ log £e + - {2eeei,r - e|) , + - 



iej 



^i) I £i 



-2 

'A 

M 

16 



+ SAej + 



l)e(,eii^r- (5.13) 



We now consider the matching of the terms with angular dependence cos {29) . The terms in 
the outer region (cf. (|5.11|) ) with such dependence are: 



eiEi^r cos (29) . 



(5.14) 



This term must be matched with the contributions due to C/4^2,2- Using (jS.Slah we obtain that we 
need to match (j5.14p with: 



1^^^ \Y\'^'-' + V2C2K, (4 



^i^i,- 



cos (261) 



The matching of (|5.14p and (jS.lSp requires: 

Ba2 



0[{e,f^+^\ asr 



00. 



(5.15) 



(5.16) 



Computing higher order terms in the outer expansion it would be possible to derive more precise 
formulas for -64,2. Basically this would require to compute higher order asymptotics of the function 

il (y) as y — )• ita. The next order correction to il in (j4.10ap is of order C\Y\ ~ cos {29) for some 
C € M. This would give exactly the behavior ()5.16p . However, since the detailed form of these 
terms will not play any role in the following, we will not continue with this analysis. The matching 
of (j5.14p and ()5.15p requires also: \/2C2i^2 = 1 and this is just a consequence of (|3.79p . 

We now consider the matching of the terms with dependence cos (40) . The term with this 
angular dependence in (j5.1ip is: 

I cos (40) 



This term must be matched with the contributions due to C/4,2,3. Due to p.l06ap the inner contri- 
bution to be matched is: 



16/^4 C3 (00) 1^^,2^5-2 



24ci (00) + V^C^Ki (^2,3)^ 



cos (40) 



Arguing as in the derivation of ()5.16p we observe that 03(00) = 0{e^ "^ ), showing that these 
terms are very small in the inner region. Taking into account (j3.79p we have: 



ci (00) 



28-3' 



(5.17) 
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which concludes the matching to this order of the functions $. We can also obtain matchings for 
the functions V. We are just interested in the first term on the right-hand side of (I3.106b|) since it 
gives a term of order one for |y| of order one. The remaining terms give contributions of order ef 
and we will ignore them. The term to be matched for V is 3(ci(oo)/e^) |y| cos (40). The gradient 
of this term with respect to y yields: 

and using polar coordinates, as well as (j5.17p this becomes: 



|yp 



26 
On the other hand, the term in (|5.ip containing cubic terms is: 



(cos (3(9), -sin (36')). (5.18) 



(a-y)|y|2a 2(a-y)'^a \y\''y (a-y)^y 

i7i6 ' i7|8 '" A I |4 TTe ' 

a a 4 a a 



which can be transformed, using polar coordinates in: 

|y|3 

J-^ (4cos^6'-3cos6',sin6l-4cos2 6'sin6l) . (5.19) 



Standard trigonometric formulas show that (j5.18p and (I5.19P are the same. 

5.6 Analysis of the ODE (15.131) and the derivation of the final profile. 

Neglecting terms of order 0{{ei^r)'^) that will be seen to have a size of order 0((e^) /r) as r — t- oo 
we obtain: 

egei^r log Si + MeiEi^r = Lej 

with 

5 3 

M=- + 8B and L= 8A. (5.20) 

4 32 ^ ^ 

Integrating this equation, we obtain: 

^(loge,)2 + 2M^ (loge,) = 2L + O ((e^,.)') , 
whence: 



(loge^)^ + 2M(loge^) = 2Lt + O f (e^,r) ) as r ^ oo, 
where we have used that log ee is of order y/r to the leading order. Therefore: 



loge£ = — v2Lr — M + o(l) as r — )■ oo. 

Then: 

Q = /3e-"v^- (l + o(l)) asr^oo, (5.21) 

where 

a^V2L = J^ - 16A, /3 = e~'' = e-^l^~^^ . 
V 16 
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In the original variable, the leading order corresponding to (j5.2ip is: 



I3^/T~i e~"Vli°g(^-*)l. (5.22) 

Notice that since ^ < 0, which we have checked numerically as was already mentioned in 
Remark 14.21 the constant L in ()5.20p is positive and then a is a real positive number. 

The asymptotics (j5.22p characterizes the width of the peaks where the mass of u is concentrated. 
The characteristic distance between these peaks is of order: 

D = AVT - t. (5.23) 

Remark 5.1 It is interesting to notice that the formulas (|5.23p provide information about the 
characteristic distance to which two peaks, with masses close to Svr, and concentrated in a width of 
order w, must be, in order to obtain blow-up with two peaks aggregating together. Notice that, for 
w small we have the following approximation for the critical distance required to have simultaneous 
blow-up and aggregation of the two peaks 



4e w ( \ 

D = — exp ( ay 2 |log (u;)| j as lu — )• 0. 



/5 

By critical distance we understand the distance at which two peaks containing a mass close to 
Svr in an area with radius w, should be localized in order to obtain singularity formation with an 
aggregating mass 16tt. 

The numerical factor Ae~"' //3 cannot be expected to be really accurate if the concentrating masses 
in the initial peaks are not distributed exactly according to the stationary solutions (12. Sp . 



Remark 5.2 Assuming that the asymptotics for (y) stated in Remark \4 ■ 5\ holds, we can obtain an 
asymptotic formula for u (x, T) as x ^ xq using the methods in ( U9^ )- Indeed, using Remark \4.5\ as 
well as (|2.1ap . (|4.1ap we can approximate u {x, t) for any t < T, i ^ T and \x — xq\ = L\/T — t, L 
large. In such regions u is basically constant in domains with a "parabolic size" \JT — i. Therefore 
the equation (II. lap can be approximated as an ODE for times i < t < T. This allows to approximate 
u (x,T) as: 



^' _._._/ ._ /,_^.^ _„,2 



U (x, r) ~ 2 6Xp —'^Ci\l log \x — Xo\ ] if (6) aS X ^ Xq. 

|x-xo| V 



It is interesting to notice that the function (p (6) mentioned in Remark \4-5\ gives the angular 
dependence of u at the blow-up point. Therefore, a more detailed study of the asymptotics of the 
solutions of (14. 4p as \y\ — t- oo would be in order. 

6 Geometric configurations of singular self-similar solutions. 

In most of the previous computations we have assumed that $ (y, r) approaches one very specific 
singular solution of ()2.3ap . (I2.3bp with the form (12. 9p . However, there exist many other solutions 
of the system (j2.3ap . ()2.3bp that could be taken as possible limits of ^ {y, r). The problem (j2.3ap . 
(|2.3bp is meaningless if we assume that <1> is just a measure, or even a sum of Dirac masses. 
However, having in mind the matching arguments in the previous sections, it is natural to assume 
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that $ has the form ()2.4p (i.e. all the masses of the peaks are Svr) and also that the equation must 
be understood as (12. 5p or, in an equivalent way, that a given peak does not interact with itself, 
something that can be justified "a posteriori" due to the local symmetry of the peaks during the 
process of aggregation. 

In this section we just obtain a few examples of solutions of (j2.5p . It is important to remark 
that the existence of these solutions does not guarantee the existence of solutions of the original 
problem ()l.lap - ()l.lbp . Indeed, although the formal arguments described in the previous Sections 
can be extended without much difficulty to more general self-similar solutions a crucial condition 
that must be satisfied, in order to obtain a meaningful equation for the width of the peaks e^, is 
the inequality: 16~^ + 2~^ — SA^ > with Ag would be a constant defined in a manner analogous 
to Lemma 14.11 for the corresponding elliptic problem. 

We do not attempt to derive a complete classification of all the solutions of (j2.5p . However, we 
will describe some particular classes of these solutions in order to illustrate the type of geometries 
that can arise during the aggregation of multiple peaks. The cases under consideration will be the 
following ones: points in a line, regular polygons, several polygons with different sizes combined, 
complete classification of solutions for A^ = 2, 3, and particular results for N = 4,5. 

We remark that the sum of the right hand side of (j2.5p vanishes for any N > 2 and for any 
configuration of points {yj} as it can be seen by symmetrization: 

TV 

Y.yj = o. (6.1) 

6.1 Solutions where all the peaks are in a line. 

We begin with solutions of ()2.5p where all the points {yj} are placed in a line. We can assume that 
this line is the horizontal coordinate axis. Then yj = {xj,0) for some real numbers {xj}:J^-|^. Then 



()2.5p becomes: 

f-4 E ,r'"I|2 =Q' ^■ = l'2,...,iV, N>2. (6.2) 



N 



=1, ^T^i 



I ^j ~ ^t- 1 



Proposition 6.1 For every integer N >2 there exists a unique solution of (j6.2p . The solution is 
invariant, up to the rearrangement of indexes, by the transformation Xj i— )■ —Xj. 

Proof. This problem can be reformulated in a variational form because the solutions of (|6.2p 
can be obtained as the minimizers of: 

N , .2 N N 

E{xi,X2,...,xn) = ^^ 2^ ^ log\xk-Xi\. (6.3) 

fc=i 1=1 k=i,e=/=k 

The functional E {xi,X2, ■■■,xj\[) is strictly convex and lower bounded in the convex set {— oo < 
xi < X2 < •.. < xn < oo}. Therefore there exists a unique minimizer where (j6.2p holds. Moreover, 
symmetry considerations prove the invariance mentioned in the statement. ■ 



Remark 6.2 The solutions of (12. Sp can be characterized in general by means of the extremal 
points of a functional similar to the one in (p.'6\) if the points {yj} are not aligned. However, in 
such general cases, the convexity properties of the functional are not satisfied and therefore, the 
functional does not allow to obtain information about the solutions in an easy manner. 
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6.2 Regular polygons. 

Proposition 6.3 For every integer N > 2 there exists a solution of ()2.5p with the points {y^} 
placed at the vertices of a regular N -sided polygon centered at the origin. The solution is unique 
up to rotation of coordinates. Moreover, the points lie on the circle with radius 2\/N — 1 centered 
at the origin. 

Proof. It is convenient to reformulate ()2.5p using complex variables. Let us write yj = 
{yj,R, yj,i) and Zj = yj^R + iyjj £ C. Then (|2.5|) becomes: 

N 

1 = 4 E r^' ^■ = ^'-^' 

or equivalently, 

N 1 

% = 8 Y. —^^ 3 = l,-N. (6.4) 



We now look for solutions with the form: 



2lTJ ■ 

z,=pe-ir\ j = l,...N, p>0. (6.5) 



Plugging (j6.5p into (|6.4p we obtain: 

■2 ^N-U l + (-l)^ , , 

+ ^^, (6.6) 

where [x] stands for the largest integer not greater than 2: € M. This equation determines p for 
each value of A^. We actually have: 

p = 2VN -1. (6.7) 



This shows that there exists a solution of (12. Sp constructing a regular A^-sided polygon. The center 
of the polygon is necessarily at the origin because of (j6.ip . ■ 



6.3 Classification of solutions for the cases N = 2 and A^ = 3. 

In these particular cases we can characterize uniquely all the solutions of (j2.5p . The problem 
becomes more complicated if the number A^ increases, because, as it will be seen later, the number 
of geometrical configurations increases with A". 

6.3.1 The case A^ = 2. 

Proposition 6.4 Suppose that N = 2. Then a solution of ()2.5p is uniquely given by yi = (—2,0), 
2/2 = (2, 0) up to rotation of coordinates. 

Proof. Due to (|6.ip we have 2/2 = — 2/1- We can assume, up to rotation, that yi = (xi,0) with 

xi > 0. Then (12. 5p is reduced to: 

xi 2 

y ~ ^' 

whence xi = 2. This simultaneously proves the uniqueness of the obtained solution in the class of 
solutions studied in Subsections 16.11 and 16.21 when N = 2. ■ 
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6.3.2 The case iV = 3. 



This case is still sufficiently simple to obtain a complete classification of the solutions. There are 
just two solutions of (j2.5p up to rotation. Either the three points are in a line as in Subsection 16.11 
or in an equilateral triangle as in Subsection 16.21 



Proposition 6.5 Suppose that N = 3. Then for every solution of (j2.5p the points {yi,y2,y3} are 
placed, up to rotation, either at the ends and the intermediate point of a segment with length being 
4\/3 or at the vertices of the regular polygon with the length of the sides being 2\/6. 

Proof. Suppose first that the three points are in a line, i.e, yj = {xj,0) for some Xj G M. 
Then the line crosses the origin due to (|6.ip and, up to rotation, the resulting solution is the 
one described by means of the minimizers of the functional E in (|6.3p . In this case, they can be 
computed explicitly. Indeed, the invariance of the solution under the transformation Xj — )• —Xj 
implies that, under the assumption xi < 3:2 < X3, we have X2 
to: 

— + ^ 

.X3 



0, xi = —X3. Then (16. 2p reduces 



^ = 4 



(2X3 



_6_ 

x-i 



whence 3:3 = —xi = 2v3. Suppose now that the three points {yj} are not in aline. We will prove 
that in this case the three points are placed at the vertices of an equilateral triangle. It is convenient 
to use the complex notation of Subsection 16.21 We may assume, without loss of generality, that 
Z3 = Z3. On the other hand, using also ()6.ip we then observe that (j2.5p becomes: 



Zl 



2Azi 



{Zl - Z2){ZI - 2:3)' 



Z2 



24z2 



{Z2 - Zi){z2 - Z^) 



(6.8) 



Due to (j6.ip we have zi ■ Z2 / 0, since otherwise the three points would be aligned against the 
assumption. Taking the absolute value of (j6.8p . we then have: 



\Zl - Z2\ \Z\ 



z?,\ 



\Z2 - Zx\ \Z2 - Z'i\ 



24. 



(6.9) 



Therefore |zi — Z3I = |2;2 — ^sl =: o" > 0. On the other hand, there is nothing special about the 
point 23 and, using the rotational invariance of (j2.5p we may replace 23 by zi and prove in a similar 
way that |z3 — zi| = \z2 — zi\ = a. Therefore the three points are at the vertices of an equilateral 
triangle and the obtained solution is the corresponding one considered in Subsection 16.21 The 
precise size of the triangle can be computed using (j6.9p as o" = 2\/6. ■ 



6.4 The case iV = 4. 

We have not obtained a complete classification of the solutions of (j2.5p if A^ = 4 but we have some 
partial results suggesting that there exist at least three solutions (up to rotation). 

Notice first that we can obtain two solutions as in Subsections 16.11 and 16.21 Actually they can 
be computed explicitly. In the case of solutions with the four peaks in a line we write: 



Xl 



-R, X2 = —9R, X3 = 9R, X4 = R 



with R > and < 9 < 1. The equation (|6.2p then becomes: 



r:' 



1 



+ 



1 



1 + 



+ 



1 



9R' 



1 1 

^29^TT' 



(6.10) 
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whence, eliminating R, we obtain after some computations: 8^^ = (1 — 0^)^, whence: 



/5-2\/6. 
since 0^ G (0, 1). Using then the first equation in (|6.1Up we obtain: 



R = 2\jV6 + 3 

and this concludes the characterization of the solution with A^ = 4 and all the peaks aligned. 

If A^ = 4 we can obtain a solution with all the peaks at the vertices of a square as indicated in 
Subsection 16.21 Using (j6.5p and (j6.6p we obtain that the vertices are at the points: 

Zj = 2V3e^\ j = 0,1,2,3. 

We remark that it is possible to obtain another solution in the case A^ = 4 that is neither of 
the ones in Subsections 16.11 nor [ 



Proposition 6.6 Suppose that N = 4. Then there exist a solution of ()2.5p with one peak at the 
origin and three remaining peaks at the vertices of an equilateral triangle. 

Proof. We look for a solution with the form: yj = pe^'^^'^''^'^, j = 1,2,3, y^ = 0. Due to the 
symmetry of solutions under the rotation of an angle 27r/3, the equation ()2.5p becomes: 

2 1 —i 1 -—i 

P'^ 1 — 63* 1 — 63' 



^ (yi + ?/2 + ys) = 0, — = 1 + 



1 — ) 

1-63' 



+ 



2tt , 
1-6"^' 



The first equation is automatically satisfied by ()6.ip . whereas the second one gives p = 4. ■ 

6.5 A^ = 5 case 

In this case we do not attempt to obtain a complete classification of the solutions, but indicate 
some examples to illustrate what type of solution can arise. We can obtain solutions with all the 
peaks in a line as in Subsection 16.11 In this case we have, due to the symmetry of the problem: 

yi = -R, 2/2 = -OR , 2/3 = 0, yi = eR , y5 = R, 
where R > and 6 G (0, 1). The equations ()2.5p are then reduced to: 

i?2 _ 3 2 OR'^ 3 26 



2 1-02' 8 26 1-02' 

Eliminating R^ we obtain 30^ — 140^ _|_ 3 — q^ whence: 



7 2^ 

3 3 



Therefore: 



R = -^/3^/W^JVT0 - 2 ('^/lO + 2 
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We can obtain also a solution where the peaks are placed at the vertices of a regular pentagon. 
Using ()6.5p and (I6.6p we obtain: 

Zj = Ae^\ j = 0,1, 2, 3, 4. 

in the complex notation. 

There is also one solution that consists of one peak at the origin and the other four peaks 
at the vertices of one square centered at the origin. Assuming that the peaks are at the points 
Zj = /3e"2"*, j = 0, 1, 2, 3, we obtain: 



Zn 



2V^e^\ j = 0,1,2,3. 



solving the equations (12. 5p . 

We finally remark that in the case A^ = 5 it is possible to obtain one distribution of peaks whose 
only symmetry is the reflection with respect to a line. More precisely, we have: 

Proposition 6.7 There exists a solution of (j2.5p with the points {yu} placed, in terms of the 
complex notation of Subsection \6.^ at the following positions: 

-Zfc = a;fc G M for A; = 1, 2, 3, Z4 = a + i/3, z^ = a — ip (6-11) 

with a < 0, /3 > 0. 

Proof. We prove the existence of a solution of (|2.5p with the form (j6.1ip by means of a 
topological argument. Due to ()6.ip we have: 

X1 + X2 + X3 
a = . (6.12) 

We assume that a is chosen as in (I6.12p . On the other hand, we can obtain an equation for f3 using 
the vertical component (or imaginary part in complex notation) of (j2.5p with j = 4: 

1 '^ 1 1 

with a given by (j6.12p . Since the right-hand side of (|6.13p is a decreasing function of /?, we see that 
there exists a unique solution of (I6.13|) with /3 > for any {xi,X2,X3) in the set —00 < xi < X2 < 
3:3 < 00. We denote it as f3 (xi,X2,X3). Moreover, notice that ()6.13p implies 

P{xi,X2,X3)>2. (6.14) 

Equations (|2.5p with j = 1, 2, 3 is reduced, due to (j6.1ip . to: 

Xk \-^ xu — Xj 2 ixk — a) , , ^ 

Y= E H 12 + , ,2^2 ' ^ = 1'2,3, (6.15) 

where a as in (I6.12p . In order to prove that there exist solutions of (IG.lSp in the cone 

C = {x = (xi, X2, X3) : —00 < xi < X2 < X3 < 00} , 
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we treat (jG.lSp as a perturbation of the equation: 

Fk{x) = ^- Y. ^^J% =0, fc = l,2,3, 

using topological degree. Since the function F [x) = {Fi, ^2,^3) (x) becomes singular at the bound- 
ary of the cone C we construct a subset U with the property that |-F(x)| > 100 on the boundary 
dU. The functions Gk (x) = 2 {xk — a) /((a — Xk) + /3^), k = 1,2, 3, are bounded in C by 2 as 
it can be easily checked considering separately the cases |xfc — a| > 1 and |xfc — a| < 1 and using 
()6.14p . Therefore we would have |G(a;)| < |i^(2;)| on dU. On the other hand, there is a unique 
nondegenerate solution of the equation F (x) = in C due to the results in Subsection 16. 11 Classical 
degree theory then shows that there exists at least one solution of {F + G) (x) = in ^, whence 
the existence of the desired solution of (j6.15p follow. 
We shall construct the subset U of the form: 

U = {x £ C : xi + e < X2, X2 + £ < X3, —R < Xk < R, k = 1,2,3} , 

where e > and R > are constants to be determined. Notice that the boundary dU is contained 
in the planes ni,2 = {x2 - xi = e} , 112,3 = {2:3 - X2 = e} , II.r = {xi = -R} , Ur = {x^ = R} . 
We will assume that 1/e is much larger than R. Along the part of the boundary dVl contained in 
the planes Ili^2, 112,3 we then have: 

^ . ^ 1 R 1 
Fi(x)> > — . 

^ ^ - e 8 - 2e 

We then proceed to consider the part of dU contained in Ur. Suppose first that x^ — xi < 1. 
Then xi > R — 1 and we obtain: 

„ , , xi 1 1 R-1 

Fi {x) = ^ + + > ^— , 

8 3:2 — xi X3 — xi 8 

which can be made larger than 100 assuming that R > 801. Suppose now that X3 — xi > 1. We 
distinguish two cases. Suppose firstly that X3 — X2 > 1. Then: 

j^ , X a;3 1 1 R 

8 X3 — xi X3 — X2 lb 

if R is large, because the last two terms are bounded by one. Suppose secondly that X3 — X2 < 1. 
Let us assume firstly that X2 — xi < 1/4. Then X3 — X2 > 3/4 and we obtain again F3 (x) > R/16. 
Suppose secondly that X2 — xi > 1/4. Then: 

^ , , X2 1 1 X2 1 ^ X2 , 

F2 ix) = ^ + > ^ ^ ^ - 4- 

8 X2 — Xi X3 — X2 8 X2 — Xi 8 

Since X3 — X2 < 1 we obtain X2 > -R — 1 and therefore F2 (x) > R/16. We then have \F (x)| > 100 
for X G dU n Ufj. The case of x G dU n U^r is similar. 

We shall observe the existence of the desired solutions of (|6.15p . It only remains to prove that 
the equation (j2.5p with j = 4 holds. This equation is just: 

3 
a sr^ OL — Xfc 

8"^(«-Xfcf + /32" 
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In order to check that this equation holds, we just notice that this is equivalent to: 

3 

a:i + X2 + ^3 _ Y^ {a- Xk) 



16 



k: 



^ (a-Xfc) +/32 



due to (|6.12|) . According to (|6.15p this equation is equivalent to: 

The last identity is trivially satisfied by syninietrization. ■ 

Remark 6.8 We have made some com,putations suggesting that in the case N = A the only trape- 
zoidal solution is the square. The only rhombic solution is also the square. Increasing the value 
of N it becomes possible to show that there are also solutions with nested squares, triangles, etc. 
However, we will not continue this discussion here. It would be interesting to determine the smallest 
number N yielding solutions without any symmetry group. 

7 Bounded domains. 

Solving the Keller-Segel model in the half circle, it is possible to obtain a wealth of shapes yielding 
aggregation at the boundary. The mass is, in all the cases ATrm with positive integers m. It 
is possible to obtain for instance Svr instead of 47r, just keeping one point at the interior of the 
domain. Notice that one must choose symmetric point configurations in order to ensure that the 
homogeneous Neumann boundary conditions are satisfied. 

Acknowledgements. The authors thank C.Cuesta, M. Fontelos, and B. Gamboa for their help 
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